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Abstract 

In this paper, based on some mathematical results obtained by Yamabe, 
Osgood, Phillips and Sarnak, we demonstrate that in dimensions three and 
higher the famous Ginzburg-Landau equations used in theory of phase tran- 
sitions can be obtained (without any approximations) by minimization of the 
Riemannian-type Hilbert-Einstein action functional for pure gravity in the pres- 
ence of cosmoloigical term. We use this observation in order to bring to comple- 
tion the work by Lifshitz (done in 1941) on group-theoretical refinements of the 
Landau theory of phase transitions. In addition, this observation allows us to 
develop a systematic extension to higher dimensions of known string-theoretic 
path integral methods developed for calculation of observables in two dimen- 
sional conformal field theories. 



1 Introduction 

1.1 Landau theory versus other theoretical methods for 



predicting crystal structure 

An article [1] in Nature written in 1988 referred to the inability, at that time, 
to predict the crystal structure of simple crystalline solids from their chemical 
composition as "a continuing scandal" in solid state physics. Recently, Oganov 
and Glass developed a method for prediction of the most stable structures along 
with some low energy metastable states for a given compound without hints 
supplied by experimental data for this compound [2]. Their method combines 
ah initio electronic structure calculations with an evolunionary algorithm. Only 
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the chemical composition is needed as an input. The method allows prediction 
of crystal strucutures at any P-T conditions. 

In view of thise remarkable achievements, it is of interest to study to what 
extent much earlier phenominological results by Landau [3] can help in solving 
the same problem of structure prediction. Landau's results were subsequently 
improved by Lifshitz in 1941 [4] and used later by Ginzburg and Landau in 
1950 [5] in connection with their study of superconductivity, as is well known. 
The Landau theory can be used only in a rather narrow domain of temperatures 
near criticality. Being thermodynamical in nature, it contains parameters whose 
actual numerical value is to a large degree arbitrary. In principle, they can be 
deduced from experimental data. The original motivation for such a theory 
came from experimental study of order-disorder phase transitions in alloys. The 
first paragraph of Chapter 14 of Vol. 5 of the famous Landau-Lifshitz course 
in theoretical physics [6] begins with description of the concept of order (to 
be defined below) in a typical CuZn alloy. For pure Cu, the lattice is face 
centered cubic (fee), while for Zn, it is hexagonal closed packed (hep) [7]. For 
an alloy composed of these two metals, the system must "decide" what kind of 
lattice to have. The decision is made based on thermodynamical considerations. 
Specifically, it is believed that if the system has reached an equilibrium at a given 
temperature and pressure, its free energy should be minimal. Since a system as 
simple as CuZn exhibits a very complicated phase diagram, depicted in Fig.l, 
such considerations are not too informative. 2 

Empirically, it is known that about 3/4 of all elements in nature are metals. 
When mixed, about 3/4 of these crystallize into one of the three most frequently 
encountered lattices with almost equal probability [8]. These are: fee, Fig. 2a), 
bec, Fig. 3 a), and hep, Fig. 2 b). In addition, Fig.l exhibits more complex 
lattices, e.g. 7 and S, whose group-theoretic description is given in Table 1. 



Table 1. 



Phase 


Composition(at % Zn) 


Symbol and # 


Prototype 


a 


0-38.27 


Fm3m 


225 


fee 





36.1 - 55.8 


Im3m 


229 


bee 





44.8 - 50 


PmSm 


221 


T? 


1 


57.0- 70.0 




217 


bee 


5 


72.45- 76.0 


P6 


174 


hep 


e 


78.0-88.0 


P&z/mmc 


194 


hep 


V 


97.17-100 


P63/mmc 


194 


hep 



A brief summary of crystallographic terminology is provided in Appendix 

2 A large number of such binary alloys, their phase diagrams and their symmetry descrip- 
tions are listed in Rcf.[8]. As compared to the textbook by Landau and Lifshitz [6], where 
actually only a small portion of CuZn phase diagram is discussed (e.g. fi and /?' phases in 
Fig.l). Rcf. [8] indicates that most binary alloys, including CuZn, exibit unexpectedly complex 
phase diagrams similar to that given in Fig.l. Subsection 6.2. contains additional information 
helpful for its understanding. 
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Figure 2: a) Basic fundamental domain for the fee lattice along with its primitive 
cell; b) the same for the hep lattice and its primitive cell 



A. The CuZn system was used in Rcf.[6] as a good example of the second order 
phase transition. A phase transition for which the order parameter changes 
continuously (with concentration) till it becomes zero at the transition point is of 
second order (or continuous) . In reality, however, such a transition takes place 
in a rather narrow concentration range (roughly, between 30 and 60 % of Zn in 
Cu), as can be seen in Fig.l. Nevertheless, this example was used by Landau 
to introduce the concept of an order parameter. To this purpose, using Fig.l, 
we fix the temperature so that the horizontal line will cross the phase boundary 
within the concentration range just mentioned. Then, on the l.h.s.with respect 
to this crossing we obtain the /3-phase depicted in Fig. 3a) where both Cu and 
Zn can occupy each lattice site with equal probability (disordered state). To 
the r.h.s.of this crossing, the occupation probabilities for Cu and Zn become 
different (ordered state). Following Landau [3], we introduce the nonnegative 
order parameter <p measuring the degree of such ordering. Clearly, it should 
be proportional to the difference between the occupation probability 1/2 in the 
disordered state (cither for Cu or Zn) and that for the ordered state which 
should be less than 1/2. As Fig.l and Table 1 indicate, the symmetry of the low 
temperature /3'-phase is that of two interpenetrating cubic lattices (as depicted 
in Fig. 3b)), while the symmetry of the high temperature /3-phase is that of 

the bec lattice and is higher. This fact is in accord with the general observa- 
tion [6] that the higher symmetry phase usually occurs at higher temperatures. 
Having said this, one has to take into account that the order parameter ip was 
defined without reference to temperature so that the order-disorder transition 
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Figure 3: The simplest example of superlattice formation: a bcc /3— type lattice 
a) is converted into /? — type lattice made of two interpenetrating Z 3 — type lat- 
tices. To facilitate understanding, a two dimensional square lattice Z 2 depicted 
in a') is splitting into two square sublattices in b') 



thus described is only concentration-dependent. This means that infinitesimal 
changes in concentration can cause changes in symmetry. 

To include both the temperature and concentration dependence, the notion 
of an order parameter should be generalized. This can be accomplished if, in- 
stead of fixing the temperature, we fix the concentration in the concentration 
range specified above. Then, a vertical line in Fig.l also will cross the phase 
boundary so that the higher symmetry phase (3 indeed occurs at higher tem- 
peratures. In view of this, to describe a phase diagram locally in the sense just 
described, Landau [3] proposed the following phenomenological expansion for 
the free energy functional T\}p\ valid in the vicinity of the transition tempera- 
ture T c : 

T{^}=To+aip + A^ 2 +Cip 3 + Bip 4 + ... . (1.1) 

It is expected, that for T > T c we are left only with F{ip} — JFo, while for 
T < T c we have to assume that coefficients a, A, etc. depend upon T, P and 
concentration c. Such an expansion is non analytic, however, since the deriva- 
tives of the free energy with respect to P and T are discontinuous at T = T c . 
This fact allows us to define the order of transition using standard rules of ther- 
modynamics. The system is undergoing a first order transition if the first order 
derivatives of F{(p} with respect to T and/or P are discontinuous at T c , while 
it undergoes a second order transition if the second derivatives of F{<p} with 
respect to P and/or T are discontinuous at T c . 
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To actually use the Landau theory under such circumstances, F{<p} is min- 
imized with respect to p. Next, the first derivatives of T{p} are calculated 
with respect to P and/or T . This minimum value of the order parameter ip 
is substituted back into the derivatives in order to check their continuity at 
T = T c . This requires making some additional assumptions on coefficients in 
the Landau expansion. For instance, for the second order transition, one should 
require a = in Eq.(l.l) while A = A(P,T) is expected to change sign at 
T = T c so that, by symmetry, it should become exactly zero at T c . This implies 
that for temperatures T close to T c the coefficient A is expected to behave as 
A(P,T) = a(P)(T - T c )/T c = ar. As for the coefficient C, it is also expected 
that it can be a function of P and T . Its presence is not mandatory, however. 
It depends upon the symmetry of the system. The details were worked out by 
Lifshitz in 1941 [4]. 



1.2 Refinements by E.M. Lifshitz 



Although the textbook by Landau and Lifshitz [6] contains only one para- 
graph describing Lifshitz refinements, there arc monographs, e.g. Rcfs. [9-11] 
and the references therein, elaborating on works by Lifshitz. Such elabora- 
tions are quite sophisticated, so we refer our readers to these monographs for 
details. Here we would like only to emphasize the key points essential for our 
work. Lifshitz improvements of Landau's theory were designed to explain phase 
transitions in superlattices. The f3 phase of CuZn mixture is an example of a 
super lattice. By definition, the superlattice, (or superstructure), is a lattice 
made of atoms of different kinds occupying existing lattice sites in an orderly 
fashion. When this happens, the lattice is subdivided into different sublattices, 
as depicted in Fig 2, so that different sublattices are occupied by different atoms 
in the alloy. To date, an enormous number of superlattices is known empiri- 
cally [7,8]. In principle, Lifshitz theory allows these predictions provided that 
the symmetry of the high temperature phase is assigned. Then, the symmetry 
of the low temperature phase (superlattice) can be predicted. Many practical 
examples of such calculations can be found in Ref.[9-ll]. It should be kept in 
mind, however, that such a theory can explain diagrams like that in Fig.l. only 
qualitatively. In its present form, it cannot make predictions about the exact 
location of the phase boundaries, which is also true of the method developed 
by Oganov. For this one needs a more detailed microscopic theory. Attempts 
to develop such theory were made (to our knowledge) only for specific systems, 
e.g. see Ref.[10,12]. It is hoped that the results of this paper and that by 
Oganov et al might stimulate further research in this area. 

The probabilistic interpretation of the order parameter ip suggests that it 
should be coordinate-independent. Yet, according to Lifshitz original work [4], 
it is allowed to be coordinate-dependent. Specifically, Lifshitz argues that the 
order parameter ip is proportional to the crystalline density p which can be 
presented as p = p + Sp, provided that p serves as the basis of representation 
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of the prescribed symmetry group Go of the high temperature phase, and 6p 
can be represented as 

n j 

The first summation takes place over all irreducible representations (except 
trivial) while the second summation takes place over the bases of these irre- 

(n) 

ducible representations. Clearly, <p\ is the analog of the earlier introduced ip 
while p corresponds to the basis of unit representation for Go- By design, dp 
is invariant under group G of lesser symmetry expected for the low tempera- 
ture phase. For brevity, following Lifshitz, we shall suppress the superscript n 
in what follows. To make a connection with Eq.(l.l) Lifshitz suggests writing 
<p t as — ifrf i where the functions r y i are chosen in such a way that 

^ 2 = E^ (1-3) 

i 

implying that X^7? = 1- Under such conditions, the expansion, Eq.(l.l), ac- 

i 

quires the following form 

T{p) = F Q +A{P, T)ip 2 W E C *P, T)ff\{li))+V A E B i& r )/j 4) ({7i})+- 

3 3 

(1-4) 

where /j 3 \ /j 4 ^ ... are the third, forth and higher order invariants of the group 
G made of j i s. The sums over j count all such respective invariants. This 
refinement by Lifshitz explains the group-theoretic nature of the coefficients in 
the original Landau expansion, Eq.(l.l). Thus, some of these coefficients are 
excluded from the expansion based on symmetry in addition to the thermody- 
namic stability requirements, discussed by Landau 3 . The results obtained thus 
far are still not too restrictive to determine the group G in the low temperature 
phase. Mathematically rigorous arguments leading to full determination of G 
can be found in the monograph by Lyubarsky [9] and are too long to be used 
for these introductory remarks. To get a feeling for these arguments we pro- 
vide some less rigorous (physical) arguments in spirit of the original paper by 
Lifshitz. 

To this purpose, we assume that the order parameter is weakly coordinate- 
dependent in the following sense. The representation of any of 230 space groups 
(Appendix A) acts in the vector space whose basis is made of functions of the 
form 

*k«(r) = UkQ e ikr . (1.5) 

For a fixed k, functions Uk a are invariant under translations in the direct lattice, 
r — > r + a, while the combination k • r in the exponent changes to k • r + k • a, 
so that if H is some vector of the reciprocal lattice, then the vectors k and 



3 It should be noted that in Landau's work these symmetry arguments are present already 
but without details. 
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k + H arc equivalent. This equivalence leads to the exponent cxp(za • H) being 
equal to 1. The subscript a in Eq.(1.5) numbers all functions with fixed k. 
Such a set of functions forms a basis for representation for the point group. The 
above expression for 'J'ka (r) can be simplified for superlattices according to the 
following arguments by Lifshitz. He noticed that lattices (expected to become 
superlattices) at higher temperature (in the disordered phase) should contain 
only one atom per unit cell A 4 . At lower temperatures the points in such a 
lattice in different cells become non- equivalent thus signalling the formation of 
a superlattice. Because of this, the number of functions u^ a is reduced to one 
(i.e. a = 1). These arguments allow us to replace as well all u^ a s in Eq.(1.5) 
by unity. Next, applying all rotations related to the discrete subgroup (of G) of 
point symmetries at fixed k creates an orbit for such a subgroup, called a star. 
Since the expansion, Eq.(1.2), is made for a given star (the number of elements 
in a star thus forming a basis), it is only natural to require the coefficient A(P, T) 
to be k-depenent. Suppose now that the low temperature phase corresponds to 
some k=ko- Thcrmodynamically, for such a phase to be stable, one should look 
at A(P, T) as a function of k. When considered as a function of k, the coefficient 
A(P, T) should possess a minimum determining k . Lifshitz recognized that such 
a condition is not constructive enough. In support of the necessity of such a 
condition the following chain of arguments can be used. If k would determine 
the structure, the density Sp would be periodic with periodicity induced by k . 
In view of Eq.s (1.2) and (1.5), such a requirement apparently looks completely 
satisfactory. Nevertheless, it is incomplete, since by itself, it is not sufficient 
for the lattice to be stable. For this to happen, one should take into account 
nearby structures with k's close to ko. Since the periodicity of Sp is built into the 
expansion, Eq.(1.2), through its dependence on ^ka, the only option allowing 
structures with nearby k's to be considered lies in making ip^ in Eq. (1.2) 
weakly coordinate-dependent. In this case the expression A(P,T)ip 2 should be 
replaced by a more general expression of the type 



where the kernel is some function of P, T symmetric with respect to interchange 
of its spatial arguments. By construction, the k-dependence of A(P, T) is built 
into the transformation properties of the kernel h(r, r'). Specifically, for h(r, r') 
consider an eigenvalue problem of the type 



The choice of the kernel h(r, r') is made in such a way that the function, Eq.(1.5), 
is identified with the eigenfunction (j)(r) in Eq.(1.7). Keeping in mind that for 
superlattices, Wka can be put equal to one and, taking into account the trans- 
lational invariance, the above eigenvalue problem can be rewritten explicitly 

4 It is known from solid state physics that the elementary cells displayed in Fig.s 2 and 3 
should not be confused with the unit cells. The difference is well explained in the book by 
Ziman [13] for instance. 




(1.6) 




(1.7) 
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as __ 

^^^(r-r')e ik - r '=A je ^ r , (1.8) 

I r' 

where summation over r' is made over vectors of the direct lattice, while that 
over I is made over the k vectors forming a star. From here we obtain, 

tmjXj ■ = 5>j™( R ) e-<km ' R (1.9) 
R 

where R = r — r'. For k = k , this lattice is expected to be stable. To check 
that this is the case, we use expansion, Eq.(1.2), for 5p in Eq.(1.6). To actually 
do so, we need to be more specific about the expansion given by Eq.(1.2). To 
this purpose we write 

8p(T) = '£ t <p i (T)e ik *-' (1.10) 

i 

where k = ko + re with re being an infinitcsimally small vector. Our order 
parameter function ip^r) can also be Fourier expanded in such a way that 

^(r) = £c ki0 (K)e-- r . (1.11) 
Using expansions Eq.s (1.10) and (1.11) in Eq.(1.6) produces 
i = EEEE c ^W^oM j dr j dv'e^+^h(r-v'y^'>r'. 

i K j k' 

(1.12) 

To simplify this expression wc note that if re's in the double integral would 
be zero, the integral would reduce to the result given by Eq.(1.9). Since the 
free energy at equilibrium is minimal, by selecting the smallest eigenvalue A(ko) 
from the set {Xj} we obtain the following result for A, 

i = A(k )££c ki0 (re)c kio (re), (1.13) 

i n 

where c k . o = c_ k . o . In arriving at this result we kept only the diagonal (in re) 
term and assumed (as Lifshitz did) that Ck i0 (re)'s are very weakly dependent 
on re. Expanding the exponents in the double integral in Eq.(1.12) in powers 
of re, keeping only terms of the first order in re, and taking into account the 
perturbative nature of such a calculation leads to the following result: 

J dr J dr'e- l ^+^- r h(r - r 'y(^+^^' 

= \ + ik ■ j dr j dr'(r' - r) e - jk * r /i(r' - r)e ik * r ' + 0(re 2 ) 

= \ 1 + ik- J dr J dr'r'h(r' - r){e- lk ' r V k - r ' - e - ik * r 'e ik * r } 

+0(re 2 ). (1.14) 
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While the terms of order 0(k 2 ) will be considered in the next subsection, here 
we would like to combine Eq.s (1.12), (1.14) in order to present our results in a 
physically suggestive form. Taking into account Eq.(1.9) we obtain: 

^ = EEi A ( k 0) +i ( K -^ A ( k ) Ik=k0 +0(K 2 )}c^,(K)cijK) (1.15) 
i n 

where, clearly, 

^A(k) | k=ko = E E v ' h{v ' - r){e-^ r e^ r ' - e~^ r ' e^ r } (1.16) 

r r' 

with integrals being replaced by lattice sums. Thus, in accord with Lifshitz [4,6], 
who only suggests (without derivation) that the requirement J^A(k) |k=k = is 
equivalent to the requirement of vanishing of invariants of the type {e~ lki ' r e lki ' r — 
g-iki-r giki-rj^ we jj ave j us ^ demonstrated that this is indeed the case. This 
equivalence is absolutely essential for determination of the symmetry of the low 
temperature phase. Its existence also implies the vanishing of the coefficient C 
in the Landau expansion, Eq.(l.l), [10]. These two conditions are the neces- 
sary and sufficient conditions for the second order (continuous) phase transition 
to take place. Violation of these conditions causes the order of phase transi- 
tion change from second to first. The details can be found in previously cited 
literature. From the preceding discussion, it should be clear that both options 
depend crucially on the assigned symmetry of the disordered phase which should 
be artificially introduced into the theory. 



1.3 Refinements by Ginzburg, Landau and Wilson. LGW 
free energy functional 

Ginzburg and Landau [5] applied the Landau theory for study of supercon- 
ductivity in restricted geometries. Subsequently (albeit in a different context) 
the Ginzburg-Landau (G-L) functional was reobtained by Kenneth Wilson, 
who not only recovered the G-L functional, but explained what approximations 
should be made in order to recover this functional from the microscopic (Ising) 
model. In view of his refinements, the G-L functional is frequently referred to 
as the Landau-Ginzburg- Wilson (LGW) functional [14]. It is written as 

Hf) = ^ + lJ <M C ( V ^W) 2 + V(x) + ^ (x) - 2fc(xMx)}, (1.17) 

where in the case of the Ising model the last term indirectly describes the cou- 
pling of the Ising spins to the external magnetic field. In the absence of such 
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a field the above functional describes the second order phase transition as dis- 
cussed earlier 5 . Usually, it is rather difficult to bring the coefficients of the 
Landau expansion in exact correspondence with parameters of respective lat- 
tice models [15,16]. This is especially true for the constant b. It is expected, 
however, that under any circumstance, b is only very weakly temperature and 
pressure dependent, and normally it is greater than zero. The case when both 
a and b are zero is a special case. It determines the so called tricritical point 
in the P,T parameter space. Physically, there can be only one tricritical point 
(P t , T t ) determined by the equations 

a(P t ,T t ) =0;6(P t ,T t ) = 0. (1.18) 

Below the tricritical temperature T t the parameter b becomes negative and the 
phase transition is of first order. In this case, it is necessary to introduce an extra 
term <?<^ 6 (x) into the LGW functional, Eq.(1.17), with g > 0. It is normally 
expected that the parameter g is practically constant below T c [17]. 

Eq.(1.17) contains the magnetic field h whose role we would like to discuss 
now. To this purpose, let us consider the extremum of the LGW functional near 
T c .We obtain, 

yffi = -c V 2 <P + aip + bip 3 - h = 0. (1.19) 

Consider a special case of a constant (coordinate-independent) field <p, <p(x) = 
tp . In this case, Eq.(1.19) is reduced to 

aip + hp% = h. (1.20) 

For a > and a weak magnetic field h, we obtain p Q — h/a. For a < we 
obtain solutions even for h = 0. These are ip = ±(|a| and ip — 0. The 
solution ip — is not a minimum of the free energy and, hence, should be 
discarded. Of the two other solutions, the system chooses one of them (which 
is interpreted in the literature as spontaneous symmetry breaking). They both 
have the same free energy 6 . Let now <p(x) = f Q + Sep = ip Q + r](x), so that we 
can Taylor series expand ^F{f}, thus obtaining 



S 2 



5ip(x.)5ip(x') h=a 



J?(x)77(x') + 

(1.21) 



where 



T{ip} \ v=v > = (-c V 2 +a + Sbipl) 5(x - x'). (1.22) 



5<p(x)6ip(x.') h=Q 

5 In this work we shall call the functional F{(p} — To (with h = 0) as Ginzburg-Landau 
(or G-L) and wc shall call it LGW if it is used in the exponent of the path integral, Eq. (1.27) 
below. Such distinction is needed for mathematical reasons : for the G-L functional the order 
parameter ip is stictly nonncgativc as wc had explained already in the main text. This allows 
us to relate the G-L and Yamabc functionals (Section 3) to each other. 

6 Once the choice is made, the order parameter can be considered in all subsequent calcu- 
lations as positive. 
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Clearly, we are interested only in the phase with lowest free energy. In this case 
a < and, therefore, we obtain, — \a\ + 3bipo = 2 \a\. This result implies that, 
indeed, the selected minimum is stable since the spectrum of the operator given 
above contains only nonnegative eigenvalues. This fact was used by Levanyuk 
and then by Ginzburg in order to determine the limits of validity of the Landau 
theory. Details of their arguments can be found, for example, in Refs[14,18]. 
Following these references, the Ginzburg number Gi is defined in 3 dimensions 
as 

Gi=6 2 T c 2 /a(P)c 3 , (1.23) 

where a(P) was defined earlier, after Eq.(l.l). In order to explain its physical 
meaning, we need to introduce two characteristic lengths. This is accomplished 
by the use of Eq.s (1.21) and (1.22). Let m 2 = a + 36^g, and consider the 
correlators defined by 

_ (L24) 

J £%(x)] exp(-/3%(x)]) 

where %(x)] = if d d r{c (v?W) 2 + mV(x)} and (3 = 1/T, with T being 
temperature (in the system of units in which Boltzmann's constant kg = 1). 
Calculation of such averages for temperatures T close to T c is described in any 
textbook on path integrals [18,19] and is readily accomplished in any dimension 
d due to the Gaussian nature of the respective path integrals in the numerator 
and denominator. In particular, for d = 3, the result is given by 

(r?(xMx')> c r = |x - x'| , r/£ » 1. (1.25) 

1 /2 

The above expression defines the correlation length £ via £ = (c/m 2 ) . Con- 
sider the volume average of such a correlator 

^ y >* I dV < ?y(xMx') >~ |j (1.26) 

where the characteristic length r is given by r 3 , ~ V" and can be estimated by 
equating this average with the previously obtained Landau result, ipQ = \a\ jb ~ 
m 2 /b. This produces the following estimate for ro: r^ 3 ~ m 4 /T c &. Consider now 
the ratio (ro/£) 6 ~ r _1 Gi, with r defined after Eq.(l.l). Clearly, the original 
Landau theory makes sense only if m 2 /b >> £ 2 /c?"q, that is, when fluctuations 
are negligible. This happens when \t\ >> Gi, i.e. for the range of temperatures: 
\T — T c \ » T c Gi. When approaching T c , this inequality is inevitably violated. 
In this case, fluctuation corrections should be taken into account as explained 
by Wilson, Fisher, Kadanoff and many others [14,18]. Theoretically, however, 
one can look at these results in different space dimensions. This allows one to 
introduce the upper d u and the lower di critical dimensions. For d > d u , the 
Landau-Lifshitz theory can be used with confidence since fluctuations become 
unimportant. For d < di fluctuations are so significant at any non zero tem- 
peratures that phase transitions do not occur. Exactly at d u fluctuations are 
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still important. Technically speaking, d u is the dimensionality at which LGW 
model is renormalizablc. This means that the functional integral of the type 

7 = yZ%(x)]exp{-/3S^(x)]} (1.27) 

with 5[y>(x)] = F{<p} — ^{(Pq} can be consistently calculated perturbativcly by 
expanding the exponent in the above path integral in powers of the coupling 
constant 6, which is dimensionless exactly at d u 7 ."Consistently" here means in 
such a way that the singularities occurring in perturbative calculations can be 
consistently removed, thus making the LGW model renormalizablc. Details can 
be found in the literature [15,19]. The sophisticated machinery of the renormal- 
ization group method, a by-product of such a renormalization procedure, allows 
one to calculate the critical exponents associated with various thermodynamic 
observables. 



1.4 Description of the problems to be solved and organi- 
zation of the rest of the paper 

The results just described allow us to formulate the following problems to be 
considered in the rest of our paper. 

1. As the Fig. 1 and the Table 1 indicate, symmetry of the high temperature 
(disordered) phase for CuZn alloy is fee as depicted in Fig. 2a. The existing G- 
L theory needs this information to be supplied in advance as an input in 
order to determine the cascade of symmetry-breaking changes taking place either 
by lowering temperature at fixed concentration or by changing concentration 
at fixed temperature. Can such information be obtained by the proper re- 
interpretation of the G-L theory? 

2. There is well a developed string- inspired formalism of conformal field 
theories in two dimensions. It takes into account full conformal invariance at 
criticality [20]. Can the G-L theory be modified in order to be in accord with 
the existing two dimensional formalism ? That is to say, in higher dimensions 
the LGW functional is normally used in the exponent of the path integral, e.g. 
see Eq.(1.27). Can LGW theory be modified (to account for full conformal 
invariance) to be used in the exponent of such a string-inspired path integral? 
To what extent does such a modification allow us to develop the exact higher 
dimensional analogs of two dimensional manifestly full conformal invariant 
theories? By doing so, under what conditions can we recover known higher 
dimensional perturbative results for the LGW model? 

In this work we provide affirmative answers to these problems. Based on 
this, we are able to consider the whole spectrum of physically interesting prob- 
lems whose mathematical description is based entirely on solutions to problems 
just formulated. Specifically, 

7 In the case of LGW functional d u = 4 while for the tricritical GLW functional, d u = 3 
[17]. 
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a) In Section 2 we perform scaling analysis of the G-L model. We argue that 
in dimensions 3 and higher such scaling analysis is insufficient to insure the full 
conformal invariance of the G-L functional even at criticality (i.e. for T = T c ). 
Nevertheless, we argue that such analysis, when properly (re) interpreted, is 
helpful in finding the correct form of a conformally invariant G-L functional. It 
is found in Section 3, where it is shown to coincide with the Yamabe functional 
known in mathematics literature for some time[21-25]. 

b) In section 4 we argue that with respect to conformal variations, this func- 
tional is equivalent to the Hilbert-Einstein (H-E) functional for pure gravity in 
the presence of the cosmological term. The physical meaning of this cosmolog- 
ical term in the context of the G-L theory is explained in Sections 3 and 4. Its 
presence is of importance for both statistical and high energy physics. 

c) To avoid confusion, only statistical physics applications are treated in this 
paper. In particular, to obtain better insight into the nature of the Yamabe 
functional (which can be used only in dimensions 3 and higher) , in Section 5 we 
describe in some detail its two dimensional analog. This Yamabe-like functional 
was considered in detail in the work by Osgood, Phillips and Sarnak (OPS), 
Ref.[26]. Section 5 is not merely a review of the OPS work. In it, we accom- 
plish several tasks. First, we connect the OPS results with those independently 
obtained in physics literature [27,28], where the OPS-type functional can be 
found in the exponents of string-theoretic path integrals. Such path integrals 
are used for calculation of observables (e.g. averages of the vertex operators) 
in CFT. They establish the most efficient direct link between the string and 
the CFT formalisms. Next, we argue that the OPS functional is the exact two 
dimensional analog of the Yamabe functional. The Yamabe should be used (in- 
stead of the OPS) in higher dimensions. This observation allows us to extend 
(later, in Section 6) the available two dimensional string-theoretic results to 
higher dimensions. 

d) In particular, as a precursor to these 3 dimensional calculations, we 
reinterpret thermodynamically the inequalities obtained in OPS paper. Such a 
rcintcrpretation is essential in determining the symmetry of the high tempera- 
ture phase thus allowing us to complete the work by Lifshitz. The results by 
Chowla and Selberg [29] were used in order to recalculate exactly the values of 
functional determinants for 2 dimensional lattices of different symmetry. The 
logarithms of these determinants are associated with the respective free ener- 
gies. The obtained exact spectrum of free energies in two dimensions provides 
us with useful reference for analogous 3 dimensional calculations done in the 
next section. 

e) Even though the analytical expression for such a 3 dimensional spectrum 
can also be obtained formally in closed form, as we demonstrate, in practice, 
some numerical calculations were required in Section 6. The results of these 
calculations are in complete qualitative accord with the exact two dimensional 
results obtained earlier. Some details of our calculations are presented in appen- 
dices A through D. 
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2 Scaling analysis and conformal invariance of 



the Ginzburg-Landau functional 



The conventional scaling analysis of the G-L functional routinely used in physics 
literature can be found, for example, in the book by Amit [15]. This analysis 
differs somewhat from that for the 4 model as described in the monograph by 
Itzykson and Zuber [30] . For the sake of uninterrupted reading we would like 
to provide a sketch of arguments for both cases now. 

We begin with the (f> 4 model. Let C(x) be the Lagrangian of this scalar 
field model whose action functional in d dimensions S[4>] is given by S[4>] — 
J d d xC(x). Let furthermore A be some nonnegative parameter. Then the re- 
quirement that S[<j>] is independent of A, i.e. J d d xC{x) = J d d x\ C(Xx), leads 
to the constraint 

J d d x(x • A + d)C(x) = (2.1) 

obtained by differentiation of S[cj>] with respect to A with A = 1 at the end of 
calculation. For C(x) given by 

C{x) = \ ( V 0) 2 + + (2.2) 

the change of C{x) under the infinitesimal scale transformation is given by 8 

Te = {X ' iL + d)C[x) +{d 4) 1 04 " m2<t)2 - (2 - 3) 

Comparison between Eq.s(2.1) and (2.3) implies that the action S[4>] is scale 
invariant if d = 4 and m 2 = 0. The result just obtained raises the immediate 
question: given that the massless G-L action is scale invariant for d = 4, will 
it be also conformally invariant under the same conditions? We provide the 
answer to this question in several steps. 

First, let M be some Riemannian manifold whose metric is g. Then, any 
metric g conformal to g can be written as g = e^g with / being a smooth real 
valued function on M [23]. Let A g be the Laplacian associated with metric g 9 
and, accordingly, let Ag be the Laplacian associated with metric g. Richardson 
[31] demonstrated that 

A s = e-/A s -i(^-l)/e-/A s -i^-l) e -/(A s /) + i^-l)e-/A s o/ (2.4) 



8 In arriving at this result we took into account Eq. (13-40) of Rcf. [30] along with condition 
D = ^ — 1, with D defined in Eq. (13-38) of the same reference. Also, we had changed signs 
(as compared to the original sourse) in C(x) to be in accord with accepted conventions for the 
G-L functional. 

9 That is Ag* = - (det g) ~ 2 di (g ij (det g) 5 dj *) for some scalar function *(rr). 
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(here (A g o /) ^ should be understood as A 9 (/(x)5'(x)). Using basic facts from 
bosonic string string theory, we notice at once that only for d = 2 does one ob- 
tain the conformally invariant (string-type) functional S[X] = J M d 2 x y /g(\/ g X)- 
(VgX) (to be discussed later in Section 5). For d > 2, in view of Eq.(2.4), the 
conformal invariance of the (f> model is destroyed. In particular, this means that 
Eq.(1.19) (for h = 0) is not conformally invariant even at criticality (a = 0)! 
Thus, using conventional field-theoretic perturbational methods one encounters 
a problem already at the zeroth order (in the coupling constant G) level, unless 
d = 2. This problem formally does not occur if one requires our physical model 
to be only scale invariant. This is undesirable however in view of the fact that 
in 2 dimensions an arbitrary conformal transformation of the metric tensor g 
of the underlying two dimensional manifold M is permissible at criticality [20] . 
Abandoning the requirement of general conformal invariance in 2 dimensions in 
favor of scale invariance in higher dimensions would destroy all known string- 
theoretic methods of obtaining exact results in two dimensions. Our general 
understanding of critical phenomena depends crucially on our ability to solve 
two dimensional models exactly. All critical properties for the same type of 
models in higher dimensions are expected to hold even in the absence of exact 
solvability. Fortunately, the situation can be considerably improved by reana- 
lyzing and properly reinterpreting the scaling results for dimensions higher than 
2. 

This observation leads to the next step in our arguments. Using the book 
by Amit [15], we consider first scaling of the non interacting (free) G-L theory 
whose action functional is given by 

S[cj>} = J d d x{(x/cpf +m 2 4> 2 }. (2.5) 

Suppose now that, upon rescaling, the field </> transforms according to the rule 10 : 

~(f>{Lx) = L"(j)(x). (2.6) 

If we require 

J d d x{{sj4>? + m 2 <f> 2 } = j d d xL d {(^~4>) 2 + m 2 A (2.7) 
and use Eq.(2.6) we obtain, 

S[<j>] = J d d xL d {(^(f ) ) 2 L 2 ^- 2 + rh 2 2 L 2 "}. (2.8) 

In order for the functional S[4>] to be scale invariant the mass m 2 should scale as: 
m 2 = m 2 L~ 2 . With this requirement the exponent lu is found to be u = 1 — |. 
This scaling of the mass is in accord with Amit [15], page 26, Eq.(2.72). It 
comes directly from the rescaling of the G-L functional, Eq.(1.17), causing the 



10 Notice that this is already a special kind of conformal transformation 



16 



coefficient of the gradient term to become one. In the standard field-theoretic 
</> 4 -type model, the coefficient of the gradient term is one already by the existing 
convention and, hence, the free field model (g = 0) is scale-invariant only if the 
mass term is zero in accord with Eq.(2.3). Hence, although the field-theoretic 
methods are applicable to both S[<f>] and the G-L functionals, the G-L functional 
is not entirely equivalent to the standard (ft 4 model functional with respect to 
scaling properties. 

The above scaling can be done a bit differently. This is also discussed in 
the book by Amit. We would like to use such different scaling to our advantage. 
Following Amit [15], we notice that although the action S[(ft] is scale invariant, 
there is some freedom of choice for the dimensionality of the field (ft. For instance, 
instead of S[(p] we consider 

S[4>] = ±i J 'd d x{( V (ft) 2 +m 2 (ft 2 } (2.9) 

with a d being some volume, say, a d = J d d x. Then, by repeating arguments 
associated with Eq.(2.8) we obtain u> = 1 (instead of previously obtained lo = 
1 — |), e.g. see Eq.(2-66) in the book by Amit [15]. Although, from the point 
of view of scaling analysis both results are actually equivalent, they become 
quite different if we want to extend such scaling analysis by considering general 
conformal transformations. Although, in view of Eq.(2.4), such a task seems 
impossible to accomplish, fortunately, this is not true, as we demonstrate below. 

The next step can be made by noticing that the mass term scales as scalar 
curvature R for some Ricmannian manifold M, i.e. the scaling rh 2 = m 2 L~ 2 
is exactly the same as the scaling of R given by 

R = L~ 2 R. (2.10) 

This result can be found, for example, in the book by Wald, Ref. [32], Eq.(D.9), 
page 446. In general, the scalar curvature R(g) changes under the conformal 
transformation g — e 2 * g according to the rule [23] 

R(~g) = e- 2 f{R(g) - 2(d - l)A g f - (d - l)(d - 2) |v 9 /| 2 } (2.11) 

where A g f is the Laplacian of / and sj g f is the covariant derivative defined 
with respect to metric g. From here we see that, indeed, for constant / 's the 
scaling takes place in accord with Eq.(2.10). Now, however, we can do more. 

Following Lee and Parker [23] , we would like to simplify the above expression 
for R. To this purpose we introduce a substitution: e 2 ^ = tp p ~ 2 , where p — -r^, 
so that g = (fi p ~ 2 g. With such a substitution, Eq.(2.11) acquires the following 
form: 

R(g) = <p 1 -*{a& g <p + R( g )<p), (2.12) 

with a = 4^5^- Clearly, such an expression makes sense only for d > 3 and 
breaks down for d = 2. But we know already the action S[X] which is both 
scale and conformally invariant in d = 2. It is given after Eq.(2.4) and will be 
discussed further in Section 5. The results of this section will allow us to obtain 
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similar actions, which are both scale and conformally invariant in dimensions 3 
and higher. This is described in the next section. 



3 Ginzburg-Landau functional and the Yamabe 
problem 

We begin with the following observation. Let R(g) in Eq.(2.11) be some constant 
(that this is indeed the case will be demonstrate shortly below). Then Eq.(2.12) 
acquires the following form 

aA g (p + R(g)tp = Rig)^- 1 . (3.1) 

By noticing that p—1 = ^| we obtain at once: p—1 = 3 (fore? = 4) and p—1 = 5 
(for d = 3). These are familiar Ginzburg-Landau values for the exponents of 
interaction terms for critical, Eq.(1.19), and tricritical G-L theories 11 . Once we 
recognize these facts, the action functional producing the G-L-type Eq.(3.1) can 
be readily constructed. For this purpose it is sufficient to rewrite Eq.(2.9) in a 
manifestly covariant form. We obtain, 

S[<p] = -7 I d d x^g{a( Vg <pf + R(9)V 2 } = (3-2) 

Minimization of this functional produces the following Euler-Lagrange equation 

aAgip + R{g)ip — Xip p ~ 1 = (3.3) 
with constant A denoting the extremum value for the ratio: 

A = mtttL = inf{/S'[y>] : g conformal to g}. (3.4) 

II r II p 

In accord with Landau theory [3], it is expected that the conformal factor ip is 
a smooth nonnegative function on M achieving its extremum value <p . Com- 
parison between Eq.s(3.1) and (3.3) implies that A = R(g) as required. These 
results belong to Yamabe, who obtained the Euler-Lagrange G-L-type Eq.(3.3) 
upon minimization of the functional S[<p] without knowlege of Landau theory. 
The constant A is known in literature as the Yamabe invariant [23,33]. Its 
value is an invariant of the conformal class (M, g) . The Yamabe problem lies in 
finding a compact Ricmannian manifold (M,g) of dimension n > 3 whose metric 
is conformal to metric g producing constant scalar curvature. 

Subsequent developments, e.g. that given in Rcf. [34,35] extended this 
problem to manifolds with boundaries and to non compact manifolds. It is 

n E.g. read the discussion after Eq.(1.18) and take into account that in the tricritical case 
d u = 3. 
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not too difficult to prove that the (Yamabe-Ginzburg-Landau-like) functional is 
manifestly conformally invariant. To this purpose, we need to rewrite Eq.(2.12) 
in the following equivalent form 

ipfR(g) = (aipA g ip + R(g)ip 2 ). (3.5) 

It can be used in order to rewrite E[(p] as follows: E\ip\ = J d d Xy/gR(g). Next, 
by noting that / d d Xy/g = J d d x y /gip p , we can rewrite the Yamabe functional 
in the Hilbert-Einstein form 

Jd d xVgR(g) 

S[<P\ = ~r (3-6) 

(/ d d xs/g) - 

where both the numerator and the denominator are invariant with respect to 
changes in the metric. This becomes especially clear if we recall that R(g) is a 
constant. 

In order to use these results in statistical mechanics, we require that the 
extremum of the Yamabe functional S[ip] is realized for manifolds M whose 
scalar curvature R(g) in Eq.(3.3) is also constant. In view of the relation 
/ d d X\fg — J d d Xy/gip p 7 it is clear that for the fixed background metric g Eq.(3.3) 
can be obtained alternatively using the following variational functional 

S[tp\ = J d d x^g{a( Vg ip) 2 + R(g)f 2 } - A J d d x^gyP (3.7) 

where the Lagrange multiplier A 12 is responsible for the volume constraint. This 
form of the functional S[ip] brings this higher dimensional result in accord with 
that appropriate for two dimensions (e.g. see below Section 5, Eq.(5.24)). 

Apart from the normalizing denominator, Eq.(3.6) represents the Hilbert- 
Einstein action for pure gravity in d dimensions. In the denominator, the volume 
V taken to power |, serves to make S[ip] conformally invariant [22], page 150. 

4 Ginzburg-Landau from Hilbert-Einstein 



In this section we analyze significance of the cosmological constant in the Hilbert- 
Einstein action for gravity from the point of view of the G-L model. To this 
purpose, following Dirac [36], consider the extended Hilbert-Einstein action 
functional defined for some pseudo Ricmannian manifold M of total space-time 
dimension d, without boundary: 

S c (g)= [ R^gd d x + C f d d xj-g. (4.1.) 

JM J M 

12 Actually, A = A| = 



19 



The (cosmological) constant C is determined by the following arguments. Let 
Rij be the Ricci curvature tensor, so that the Einstein space is defined as a 
solution of the following vacuum Einstein equation 

Rij = Xgtj (4.2) 

with A being constant. From this definition it follows that 

R = dX. (4.3) 

Following Dirac [36], variation of the action S c (g) produces 

Rij-\9ijR+\cgi 3 - = o. (4.4) 

Combined use of Eq.s(4.3) and (4.4) produces as well, 

C = X(d-2). (4.5) 
Again, using Eq.s(4.3) and (4.5) we can rewrite Eq.(4.4) as follows 

R t] - X -g l0 R +^(d- 2)R 9ij = 0. (4.6) 

These results remain unchanged if we use the Riemanninan manifold M in- 
stead of pscudo Riemannian. This observation allows us to use many facts from 
the Yamabe theory [22-24], originally developed for Riemannian manifolds. 



In view of this, we argue that Eq.(4.6) can be obtained as well by varying 
the Yamabe functional, Eq.(3.6). Indeed, following Aubin [22] and Schoen [24], 
let t be some small parameter labeling the family of metrics gtj (t) = gtj + thij . 
Then, these authors demonstrate that 



dR t 
dt 



vV hij-^Vjht-RVhij (4.7) 



t=o 



and 



J^^lVW^^ (4.8) 

where, as usual, \g\ = |det gij\. 

Consider now the Yamabe functional, Eq.(3.6), but this time, written for 
the family of metrics. We have 

H(g(t)) = (V(t))T I R{g(t))DV{t) (4.9) 

JM 

where the volume is given by V(t) — J M d d Xy/ g(t) and, accordingly, DV{t)=d d 
Using Eq.s(4.7) and (4.8) in Eq.(4.9) and taking into account that the combi- 
nation y 1 V J hij — Vj K is the total divergence, produces the following 
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result: 



= V(0)- 2 ^[[ (R^/2-R^ hlJ DV(0) ( DV(0) 

~(\ -hi DV(0) I W'W(O)]. (4.10) 
z a Jm Jm 



I M JM 

If the metric g is the critical point of lZ(g(t)), then 



/ W(0) + (!-!)(/ BDV(0)) 9(j =0. (4.11) 



From here, multiplication of both sides by g u and subsequent summation pro- 
duces at once 

R-^d+(^-^) <R>d = 0, (4.12) 

where < R >= pijy J RDV(0) is the average scalar curvature. Eq.(4.12) can 
be rewritten as R =< R > . But this condition is exactly equivalent to the Ein- 
stein condition, Eq(4.2), in view of Eq.(4.3)! Hence, under such circumstances, 
Eq.s (4.11) and (4.12) are equivalent. For the set of metrics of fixed conformal 
class (i.e. related to each other by g — e 2 -' g) the variational problem for the 
G-L functional, Eq.(3.7), is exactly equivalent to the variational problem for 
the Hilbcrt-Einstein functional, Eq.(4.1), for gravity field in the vacuum in the 
presence of the cosmological constant. Moreover, the result R =< R > jsut ob- 
tained coincides with the proviously obtained Eq. 92.12) (in view of Eq. (3.4)). 
This equivalence is of major importance for application discussed in Sections 5 
and 6. 



The results displayed above become trivial for d = 2. Physically, however, 
the case d = 2 is important since it is relevant to all known statistical mechanics 
exactly solvable models treatable by methods of conformal field theories. Hence, 
we would like to discuss some modifications of the above results required for 
treatment of two dimensional analogs of the G-L theory. 



5 Ginzburg-Landau-like theory in two dimen- 
sions 



5.1 General remarks 



From field-theoretic treatments of the G-L model [19] we know that straight- 
forward analysis based on asymptotic e— expansions from the critical dimension 
(4) to the target dimension (2) is impractical. At the same time, known re- 
sults for CFT and exactly solvable models arc useful thus far only in d = 2. 
The question arises: is there an analog of the G-L Eq.(3.1) in two dimensions? 
And, if such an analog does exist, what use can be made of it? In the rest of 
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this section we provide affirmative answers to these questions. We demonstrate 
that: a) indeed, a two dimensional analog of the G-L equation does exist and is 
given by the Liouville Eq.(5.f9) below, b) the fuctional, Eq.(5.14), whose mini- 
mization produces such an equation is the exact two dimensional analog of the 
G-L-Yamabe functional, Eq.(3.2), c) these results can be (re)obtained from the 
existing string-theoretic formulations of CFT developed entirely independently, 
d) in view of the noted correspondence, these string-theoretic CFT results can 
be extended to account for the Lifshitz-type problems discussed in the Introduc- 
tion. This allows us to obtain a positive answer to the 1st Problem formulated 
in the Introduction. 

5.2 Designing the two dimensional G-L-Yamabe functional 

To discuss topics related to items a) and b) just mentioned, we begin with the 
observation that in two dimensions, Eq.(2.4) acquires a very simple form 

A~ g = e~ 2f A g , (5.1) 

where we use a factor of 2 to be in accord with Eq.(2.11) for scalar curvature. 
According to Eq.(2.11), the scalar curvature in two dimensions transforms like 

R(g) = e- 2 f{R(g) - A g 2f} (5.2) 

while the area dA = d 2 x^fg transforms like 

dA = e - 2f dA. (5.3) 

These facts immediately suggest that the previously introduced action func- 
tional 

5[X] = / d 2 xV?( Vs X) • ( Vs X) (5.4) 

is conformally invariant. Using results by Polyakov [37] and noting that (y 9 X) • 
(VgX) = g a Pd a XV j dpX fl , we need to consider the following path integral 13 

cxp(-H9)) = J D[cf>] c M-\j M d 2 x s /gg a ^d a cj ) d^) (5.5) 

where T{g) is the "free energy" 14 . Fortunately, this integral was calculated 
by Polyakov [37] for two dimensional manifolds M without boundaries and by 
Alvarez [38] for manifolds with boundaries. In this work we shall be mainly 
concerned with manifolds without boundaries. Although in the original work by 

13 Without loss of generality, we would like to consider the case of a one component field <f> 
only 

14 Usually, instead of F(g) one writes T{g)lk B T,where T is the temperature and ks is the 
Bolzmann's constant. In the present case the problems we are studying do not require specific 
values for these constants. For this reason they will be suppressed. Also, one should keep in 
mind that the free energy is always defined with respect to some reference state. This will be 
the case in our calculations as well. 
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Polyakov, one can find the final result of calculation of the above path integral, 
the details of this calculation can be found only elsewhere. In particular, we 
shall follow pedagogically written papers by Weisberger [39,40] and Osgood, 
Phillips and Sarnak [26,41, 42] (OPS). 

To begin, let g a p be some reference metric, and let g a p be a metric confor- 
mally related to it, i.e. g a p = exp(—2(p)g a p. Should the above path integral be 
for the flat (i.e.g a f3=d a p) two dimensional manifold, one would have at once the 
result: T= \ lndet Ao, where the prime indicates that the zero mode is omitted. 
Because to assume such flatness in general is too restrictive, it is appropriate 
to pose a problem: how is the path integral for the metric g related to that for 
the metric gl The paper [26] by OPS provides an answer, e.g. see Eq.(1.13) of 
this reference. To connect this equation with the free energy, we replace it by 
the equivalent expression 

ln ^ _ ln («±y->£f iM + R(9M] (5 ,, 

M 

useful in applications to strings and CFT, e.g see Ref.[43], page 637. 

It is worthwhile to provide a few computational details leading to Eq.(5.6). 
Therefore, let ip i be eigenfunctions of the Laplacian A g with eigenvalues A^ 
arranged in such a way that = Ao < Ai < A2 < • • •, i.e. 

-A g fa + Xfa = 0. (5.7) 

Then, we construct the zeta function 

00 

C 9 M=EV S (5-8) 
i=i 

in such a way that 

dot 'A g = exp(-Cj(0)) (5.9) 
with C 3 (0) = ( j^C g ( s )) s=0 • Using Eq.(5.1), we obtain as well 

-e- 2 *A g iP t + AVu - 0. (5.10) 

In particular, for a constant (p = (p we obtain 

00 _ 

i=l 

Use of Eq.(5.9) in Eq.(5.11) produces: 

C' g (0) =C>) + 2^(^-1). (5.12) 
This result was obtained with help of the known fact, Ref.[26], Eq.(1.9), that 

C s (0) = ^ - 1 (5.13) 
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with x(M) being the Euler characteristic of two dimensional manifold M with- 
out boundaries. In view of the definition, Eq.(5.9), and the Gauus-Bonnet the- 
orem, we observe that Eq.(5.6) is reduced to Eq.(5.12) for the case of constant 
conformal factor ip = <p, as required. 

It should be mentioned that Eq.(5.6) was used earlier in finite size scaling 
calculations [44-46]. Unfortunately, the authors of [45,46] did not take into 
account that the above formula is valid, strictly speaking, for manifolds M 
without boundaries while the topology considered in these works was that for 
the punctured disc and/or annulus. Rigorous results for such topologies were 
obtained by Weisberger [40] and more general case of surfaces of higher genus 
with boundaries and/or punctures was discussed in detail by OPS [42]. To 
focus on the main goals of this paper, we shall not discuss the finite size scaling 
analysis further. 

Using Eq.(5.6) and, following OPS [26], we would like consider the related 
functional defined by 15 

Hf) = \j dA g {\ Vg <p\ 2 + R(g)tp} - ttx(M) In J dA g e 2 * . (5.14) 

M M 

This functional is the exact analog of the Yamabe functional, Eq.(3.2), in di- 
mensions 3 and higher as we shall demonstrate shortly below. In the meantime, 
in view of Eq.(5.6), it can be rewritten in the following equivalent form 16 

F{<p) = -67rlndet 'A g + 7r(6 — x(M)) In A. (5.15) 

Let a be some constant, then 

F(<p + a) =J 7 (v)- (5-16) 

This signifies that the above action is scale invariant. This property is in com- 
plete accord with Eq.s (2.7) and (2.8) discussed earlier. If we impose the con- 
straint (fix the gauge): A = 1, then we end up with the Liouville-like action 
used in CFT. We would like to explain this in some detail now. 

Following OPS, it is convenient to replace the constraint A = 1 by the 
alternative constraint on the field <p: 

J (pdA g = 0. (5.17) 

M 

To demonstrate that such an imposed constraint is equivalent to the requirement 
A = 1, we note that, provided that the field ip minimizes ^(tp) subject to the 



15 Whcn comparing with the OPS paper, it should be noted that OPS use the Gaussian 
curvature K while we use the scalar curvature R = 2K 

16 Here one should understand the word "equivalence" in the sense that both functionals 
produce the same critical metrics upon minimization. 
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constraint Eq.(5.17), the field 

ip = ip - i In J exp(2iP)dA g (5.18) 

M 

minimizes !F(<p) subject to the constraint A = 1. Using these facts, minimization 
of ^(tp) produces the following Liouville equation 

M 

Comparing this result with Eq.(5.2) and taking into account that 

2n X (M) 



J exp(2ip)dA g 

M 



R(g) = const, (5.20) 



we conclude that, provided that the background metric g is given so that the 
scalar curvature R(g) (not necessarily constant) can be calculated, the Liouville 
Eq.(5.19) is exactly analogous to the previously obtained Eq.(3.3), or equiva- 
lcntly, Eq.(2.12). In view of this analogy, use of Eq.s (3.4) and (3.6) as well as 
Eq.s(5.14),(5.15) causes the functional F(<p) to attain its extremum for metric 
g of constant scalar curvature R(g). To decide if the extremum is minimum or 
maximum, we have to consider separately cases x(M) > and x(M) < 0. In 
the case of x(M) > we have only to consider manifolds homeomorphic to S 2 
so that x{M) = 2. Fortunately, this case was considered in detail by Onofri 
[47] 17 . Using his work, the following inequality 

In J dA g exp(i>) < J dA- g ii> + ^ J dA g \^\ 2 (5.21) 

s 2 s 2 s 2 

attributed to Aubin [22] and inspired by previous results by Trudinger and 
Moser [48] , is helpful for deciding whether the obtained extremum is minimum 
or maximum. Here g is the metric of the unit sphere S 2 with constant Gaussian 
curvature 1. The metric g conformal to g is given by g — cxp(2tp)g , with ip 
obeying the Liouville equation (just like Eq.(5.19)), where both R(g) and R(g) 
are constant by virtue of the initial choice of g. By combining Eq.s(5.6),(5.14) 
and (5.15) with the inequality (5.21) and taking into account that by design 
\nA g = we obtain, 

_37rln de7A^ = 4 / d M\VgM 2 + ^)-^ dA 5 cxp(2ij)>0, (5.22) 

s 2 s 2 

with equality occurring only at the extremum ip = tp* , with the function tp* 
being a solution of the Liouville Eq.(5.19). It can be shown [47] that: a) such a 



^Generalization of these two dimensional results to higher dimensional manifolds of even 
dimensions can be found in the paper by Bcckner [48] . 
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solution involves only Mobius transformations of the sphere S* 2 and that, b) the 
functional T(ip) is invariant with respect to such transformations. The case of 
x{M) < is treated in Section 2.2 of the OPS paper, Ref.[26] and leads to the 
same conclusions about extremality of the functional T(ip). 

The two dimensional results just obtained are in accord by design with results 
obtained in higher dimensions (discussed in Sections 3 and 6). In particular, the 
functional T(ip) is the exact analog of the Yamabe functional S[ip], Eq.(3.2). 
Since in both cases the functionals are "translationally" (actually, scale) in- 
variant, e.g. compare Eq.(3.6) with Eq.(5.16), in both cases, the extremum is 
realized for metric conformal to the metric of constant scalar curvature, e.g. 
compare Eq.(3.3) with the Liouvillc Eq.(5.19). 



5.3 Connections with string and CFT 



The results just obtained allow us now to discuss topic c) listed at the begin- 
ing of this section. Using known facts from string and conformal field theories, 
[27,28], it is of interest to consider averages of the vertex operators 

(f[exp(/3^)))= [ D[cf>}cxp{-S L (cf>)}f[cxp(PM^)), (5-23) 

\i=l / •* i=l 

where the Liouville action Sl{4>) is given (in notation adopted from these ref- 
erences) by 

S ^) = T- I d M\V<f>\ 2 -Qii(5)^ + 87r/2exp(a+^)]. (5.24) 

The actual values and the meaning of constants Q, Ji and a + are explained 
in these references and are of no immediate use for us. Clearly, upon proper 
rescaling, we can bring Sl(4>) to the form which agrees with ZF{ip), defined by 
Eq.(5.14), especially in the trivial case when both x{M) and jj, are zero. When 
they arc not zero, the situation in the present case becomes totally analogous 
to that discussed earlier for the Yamabe functonal. In particular, in Section 
3 we noticed that the G-L Euler-Lagrange Eq.(3.3) can be obtained either by 
minimization of the Yamabe functional, Eq.(3.2) (or (4.9)), or by minimization 
of the G-L functional, Eq.(3.7), where the coupling constant A plays the role of 
the Lagrange multiplier enforcing the volume constraint. In the present case, 
variation of the Liouville action Sl(</>) will produce the Liouville equation, e.g. 
see Eq.s (5.19)-(5.20), which is the two dimensional analog of the G-L Eq.(3.3). 
This variation is premature, however, since we can reobtain f{y>) exactly using 
the path integral, Eq.(5.23). This procedure then will lead us directly to the 
Liouville Eq.(5.19). 
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To this purpose, we need to consider the path integral, Eq.(5.23), in the 
absence of sources, i.e. when all j3 i — 0. Following ideas of Ref.s [27,28], we 
take into account that: a) 

^ J dA g R(g) = X (M) = 2-2g (5.26) 

M 

with g being genus of M and, b) the field (f> can be decomposed into 4> = 4>o + V 5 m 
such a way that 4> is coordinate-independent and ip is subject to the constraint 
given by Eq.(5.17). Then, use of the identity 

oo 

/dx exp(ax) exp(— bexp(~/x)) = —b~~T( — ) (5-27) 
7 7 

— oo 

(with T(x) being Euler's gamma function) in the path integral, Eq.(5.23), re- 
quires us to evaluate the following integral 



I = 



Q - f 

i o exp(0 o — x(M))exp((^ / dA- g cxp(a + <^)) cxp(a + ^ ) 

M 



T(-s) 



M 

with s given by 



(p, / dAq ex.p(a + <p)Y (5.28) 
"+ J 



2a 



Q X(M). (5.29) 



Using this result in Eq.(5.23), we obtain the following path integral (up to a 
constant) 

Z[<p]= J DMexp(-^)) (5.30) 
with functional !F{ip) given by 

= S L (p; M = 0) - 2^"X(M) In^ f dAg exp(a +¥ >)]. (5.31) 

M 

This functional (up to rescaling of the field if) is just the same as the Yamabe- 
like functional !F(ip) given by Eq.(5.14). 

Define now the free energy T in the usual way via T = — ln.Z[<^] (as was done 
after Eq.(5.5)) and consider the saddle point approximation to the functional 
integral Z[<p]. Then, for a spherical topology, in view of Eq.(5.22), we (re)obtain 
T > with equality obtained when ip = ip* . Sources (or the vertex operators) 
can be taken into account also, especially in view of the results of [42]. 18 



18 Since the results of this reference are useful as well as well for extension of the finite size 
scaling results given in works by Affleck [46] and Blote et al [45] , the detailed analysis of this 
case is left for further study. 
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To better understand the physical significance of the obtained results, it 
is useful to rcobtain them using a somewhat different method. The results 
obtained with this alternative method are also helpful when we shall discuss 
their higher dimensional analogs in the next section. To this purpose we write 
(up to a normalization constant) 



oo 

J D[4>]exp{-S L (4>)} = J dAe^ A Z L (A) 



(5.32) 



where 



Z L (A) = 



I Dm I ■ 

J Jm 



dAg exp(a + (f>) - A) exp(S L ((f>; \x = 0)). 



(5.33) 



If, as before, we assume that <b = 



ip, then an elementary integration over 



b produces the following explicit result for Zl{A): 



Z l (A) = —A" I ' D[tp] ( dA a exp(a+p) 
a + J Um 



1-(w+l) 



where the exponent u> is given by 



X(M)Q 
2a + 



1. 



exp(-S L (<p;fj, = 0)), 

(5.34) 

(5.35) 



Finally, using Eq.(5.34) in Eq.(5.32) produces back Eq.s (5.30) and (5.31) (again, 
up to a constant factor). An overall "-" sign can be removed by proper normal- 
ization of the path integral. These results can be used for computation of the 
correlation functions of conformal field theories (CFT). Details can be found in 
[28]. 



5.4 Completion of the work by Lifshitz in 2 dimensions 

Since in this work our main interest is investigation of higher dimensional 
analogs of the results just obtained, no further computational details related 
to two dimensional CFT will be presented in this work. Instead, we would like 
to complete our investigation related to item d) mentioned at the begining of 
this section. It will enable us to develop similar treatments for higher dimensions 
to be considered in the next section. 

We beging with combining inequalities given by Eq.s(5.21) and (5.22). This 
produces, 

-lndet 'A g - (-lndet 'A s ) > 0, (5.36) 

with equality taking place only when the metric g is equal to that for the unit 
round sphere S 2 , i.e. to g. This means that 

det 'A g < det ' A- g (5.37) 
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implying that the determinant of the Laplacian for the round sphere provides 
the upper bound for determinants of Laplacians whose metric is conformally 
equivalent to that for the round sphere S 2 . Calculation of the determinant of 
the Laplacian for the round sphere can be done with help of Eq.s(5.8) and (5.9). 
Specifically, for this case we need to calculate Z'-(0) where 

This calculation is rather difficult and can be found, for example, in [49] along 
with its multidimensional generalization. For S 2 the final result is given by 

Zi(0)=<(-1)-1 (5.39) 

From here we obtain, det ' A g < cxp{i — 4£'(— 1)}. We deliberately avoid long 
discussion leading to this result since we are more interested in similar calcu- 
lations for torus topology. In this last case the reference metric g is that for 
the flat torus T 2 = C/A, where C is the complex plane and A is some lattice. 
Evidently, x(T 2 ) = 0, so that R(g) — 0. These facts produce at once 



■In det 'A s -(-lndet 'A- g ) = ^- J dA g | Vs ¥>| 2 > 

M 



(5.40) 



with equality when ip = const. 

From the point of view of applications to statistical mechanics (or better, 
to CFT), the above inequality is useful when calculating the path integral, 
Eq.(5.30), using the saddle point method. In this case, since T{p) plays a role 
of the LGW free energy, in accord with general requirements of thermodynamics, 
the free energy attains its minimum at equilibrium, i.e. for ip = ip*. Here ip* is 
a solution of the Liouville Eq.(5.19) in accord with previously obtained result, 
Eq.(5.22), for the sphere. In the present case of toral topology, it is known 
that x{M) = 0, so that ip* = const is an acceptable solution. Such a two 
dimensional result is formally in complete accord with the higher dimensional 
G-L result, Eq.(1.20), ( h = 0) of Section 1. In the present case, the choice 
of the const is dictated by Eq.(5.18) 19 . For such a choice, F(tp*) = 0. This 
is clearly the lowest possible value for the free energy, which in G-L theory 
corresponds to the free energy at criticality (i.e. at T = T c ). Expanding 3F(ip) 
around the equilibrium value of the field ip in the path integral in Eq.(5.30) 
produces quadratic and higher order terms as usual [19]. If one ignores terms 
higher than quadratic, one obtains the standard Gaussian-type path integral 
discussed in the Introduction. It plays a major role in many two dimensional 
CFT models [20]. Justification of such a truncation is discussed in the next 
section. In the meantime, we would like to discuss calculation of this type of 



19 Prom the arguments related to Eq.(5.18) it should be clear that the restriction A = 1 is 
not essential but convenient 
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path integral in some detail. Since the path integral is Gaussian, its calculation 
is essentially the same as calculation of det ' A g . In view of Eq.(5.9), such a 
calculation involves use of the zeta function of the following type 




1 



(5.41) 



(47T2 |1| 2 ) 



This result can be easily understood using some basic information from solid 
state physics [13]. Indeed, the eigenfunctions /(l) of the Laplace operator for 
T 2 are given by /(l) = exp(z27r(^ini + I2TI2)) with vector 1 = {li, I2] <G A* being 
some vector of the dual (reciprocal) lattice A*, while the numbers {ni,n 2 } are 
related to some vector of the direct lattice A. Accordingly, the corresponding 
eigenvalues are 4-7T 2 |1| 2 20 . 

The method of calculation of Z g (s) presented in Ch 10 of Ref.[20] while being 
straightforward is not too illuminating, especially if one is contemplating its 
extension to dimensions higher than 2. We follow therefore the approach taken 
by OPS, Ref.[26], where the first Kronecker limit formula is used for evaluation 
of Z'~(0). To avoid duplications, our arguments (leading to the same results) are 
somewhat different than those used in the OPS paper. These arguments allow 
us to make additional useful connections with some facts from number theory. 21 

It is well known [51] that for each torus T 2 , the modular lattice A is given 



of the periods should be complex. Different tori are related to each other by a 
modular transformation of the type 



provided that ad — be = 1 with a, 6, c and d being some integers. The require- 
ment ad — be = 1 guarantees that the inverse transformation (i.e. t' — > r) also 
looks like Eq.(5.43), with respective constants being also integers. Invariance 
of the observables of CFT with respect to modular transformations is widely 
emphasized [20]. However, such invariance is too broad and leads to some 
inconsistencies in CFT and string theories discussed in Ref.[50]. These incon- 
sistencies can be removed if we restrict values of the modular parameter r to 
those originating from complex multiplication (CM). The concept of CM can 
be easily understood based on the following example. Consider two lattices A 
and A' such that there is a matrix A 



20 Since this result can be extended immediately to 3 and higher dimensions it will be used 
in the next section as well. 

21 Some of these facts have been presented already in our recent work, Rcf.[50]. 





(5.44) 
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so that 



u 2 = au2 + bui, (5.45) 

bJ 1 = CU>2 + du>l. 

Clearly, if we require ad — be — 1 and then form a ratio of the above two 
equations, we would obtain Eq.(5.43). This time, however, we would like to 
keep the requirement ad — be — 1 but avoid forming the ratio. This move is 
motivated by the fact that by making this ratio, we would loose the option: 
lo 2 = oluj 2 and lo 1 = auj\ for some, yet unknown, a. Substitution of these 
relations into Eq.(5.45) leads to the following eigenvalue problem 

a 2 - atrA+ det A =0. (5.46) 

But, we know already that det A =1. Hence, Eq.(5.46) can be rewritten as 

a 2 - a(a + d) +1=0. (5.47) 

This is the equation for an integer in the quadratic number field. It is easy 
to demonstrate that such an integer must belong to the imaginary quadratic 
number field. To prove this fact, we write for the roots 

ai,2 = ^±^v / (« + ^) 2 -4. (5.48) 

In order for a to be an integer belonging to the imaginary quadratic field, the 
following set of options should be explored first. These are: 

a) a = d = 0, thus producing ct\.2 = ±i; 

b) a = ±1, d = (or a = 0, d = ±1), thus producing ai ; 2 = ^(±1 ± \f— 3); 

c) a = d = ±1, thus producing a\ t 2 = ±1- 

Before analyzing these results it is helpful to recall the following definition 
from CM. A torus T 2 admits CM if it admits an automorphism A — > cA so 
that cA C A for some c ^ Z. The case when cA C A requires us to replace 
Eq.(5.48) by 

a ia = ^± 1 -^{a + dy-An (5.49) 
where the nonnegative integer n should be strictly greater than one. 

To understand the physical implications of these results, recall that if a(A) 
denotes the area of the period parallelogram associated with A, then it can be 
shown [51] that 

a(A) = - \uJ1LJ2 — 0J2&i\ ■ (5.50) 

If we now rescale oj's: u 2 = uj 1 = au>i, this produces: a(aA) = |a| 2 a(A). 
If we require that the area upon rescaling remains unchanged, we are left with 
| a | 2 = 1. This is the only acceptable option in view of the constraint A = 1 on 
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the area, which was imposed earlier. Above we obtained two types of integers 
(also units) of the imaginary quadratic field : a) the Gaussian units: ±l,±i of 
the Gaussian ring Z[i] which are the 4-th roots of unity, i.e. 1, i, i 2 = —1 and 
i = —i, and, b) the units of the ring Z[j] which are the 6-th roots of unity, 
i.e. 1, — 1, j, —j, j 2 and — j 2 , where j = 5 (— 1 + V"^3). Since these are units 
in the respective rings, automatically we get \a\ 2 = 1. In view of the constraint 
ad— be = 1, no units from other imaginary quadratic fields can be used. Indeed, 
in the most general case of CM, we have to use Eq.(5.49) instead of Eq.(5.48). 
This is equivalent to replacing the constraint ad — be = 1 by ad — be = n and 
leads to a(aA) = na(A). Clearly, since by definition n^l, for physical reasons 
we are left only with the units of the two rings just discussed. 

CM along with the area constraint is sufficient for determination of all reg- 
ular lattices in two dimensions. These are: the square (associated with i), the 
hexagonal (associated with j) and the triangular (dual of the hexagonal). These 
are described in Rcfs[52, 53]. No other translationally invariant lattices exist 
in two dimensions [53] . 

Number-theoretically, these results can be reformulated as follows (e.g. see 
Appendix B). For the Gaussian lattice A, we can take as basis Aj = [1, i] while 
for the cube root of unity (hexagonal) lattice Aj = [I, j]. It should be clear that 
the usual complex numbers z = x + iy belong to Aj[53] so that the Gaussian 
integers are made of x and y being integral. Analogously, the "hexagonal" 
numbers are made of z = x + jy with " hexagonal integers" made accordingly 
from the integral x's and y's. The norm N(z) can be defined now as is done in 
complex analysis, i.e. 

N(z) = zz. (5.51) 

The norm has a very useful property which can be formulated as follows. If 
z" = z'z, then N(z") = N(z')N(z). In particular, if z is the unit of the 
complex imaginary quadratic field, then there is some z' in the same field such 
that zz' = 1. This produces N{z')N(z) = 1, so that N(z) = 1. In the case of 
a Gaussian field this leads to the equation n 2 + m 2 = 1, producing 4 Gaussian 
units: ±l,±i. At this point it is useful to keep in mind that for the field of 
real numbers there are only 2 units:±I. In view of the results just presented 
we would like to rewrite the zeta function in Eq.(5.41) in the number-theoretic 
form. We obtain, 

^)=(i^z; eA .^ ^ 

Apart from the factor (47r 2 ) s , the obtained result corresponds to the Dedekind 
zeta function. In the case of other (higher than quadratic order) number fields, 
the situation is more complicated as explained in the book by Terras [54] , so that 
other methods should be used. These will be discussed in the next section. In 
the case of two dimensions, connections between the exactly solvable statistical 
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mechanical models and number theory have been known for some time 22 All 
these connections are based on the known fact that the Dedekind zeta function 
is related to the Dirichlet L-function via 

T' — ^-=K(*)f^ (5.53) 

^zGA* N(z) s ^ y ' ^ 71 s V ^ 

where /x is the number of units in the number field and x( n ) i s the Dirichlet 
character. In particular, for the Gaussian number field we obtain 

E' T^T- = 4 C(*) E — (5-53) 
^zeA* N(z) s w ^ n s v ; 

■' v ' n=l 

with 

f if n is even 

X 4 N = 1 , !N(n-l)/2 ., . (5-54 

[_ (—1) 11 n is odd 

while for the hexagonal number field we get 

E' Tr^r = 6 CW E — (5-55) 

v ' n=l 

with 



if n = mod 3 
X 3 (n)={ lifn=lmod3 . (5.56) 
-1 if n = —1 mod 3 

To obtain 4-Zg(s) \ s —o using Eq.s (5.53)-(5.55) is possible, but is not very 
illuminating. It is much better to use the 1st Kroneckcr limit formula valid for 
any quadratic number field [55]. In order to reveal the meaning of this formula 
and its relevance to our needs a few steps are required. First, we note that the 
area form given by Eq.(5.50) can be equivalently rewritten as 

1 2 

a(A) — — \u)\0J2 — = |wi| Imr (5.57) 

where r = ^ . Second, if 1 is the vector of the reciprocal lattice A* , then it can be 
decomposed as 1 = hcoi + I2OJ2 with l\ and I2 being some integers. Accordingly, 



|2 



(hoJl + l 2 0J2)(h UJ l + I2OJ2) 



= iLu^ih+T^ih+fh). (5.58) 

Next, since we have fixed the area a(A) so that in our case, it is equal to one, 
we obtain 

Imr = (5-59) 

l w i| 



22 These can be found for example at the Google database, e.g. sec "Number theory, Physics 
and Geometry", Les Houches, March 2003. 



33 



By denoting Imr = y we can rewrite Eq.(5.41) as follows 

Zg(s)=Y' 1 -— T = —±- T -Y' (5.60) 

(4^ |1| 2 ) (2^) 2sZ ^ |Z 1+r Z 2 | 2s 

In mathematics literature [55] the derivation of the Kronecker 1st limit formula 
is given for the function 

^' s )=Ei , ,, /. ,2. ( 5 - 61 ) 
*-^h,h + T l 2 \ 

where the prime indicates that the summation takes place over all integers 
(h,h) 7^ (0,0). For many applications, it is more convenient to consider the 
combination 

E*(T,s) = n- s T(s)±E(T,s) (5.62) 

which possess the nice analytical property [56] 

E*(t,s) =E*(t,1- s). (5.63) 

Kronecker found the Laurent expansion of E*(t,s) near s = 0. It is given by 
[57] by 

E*(t,s) = ^ + 1- In 20Fy - In \ V (r) | 2 + O(s) (5.64) 
where 7 is Euler's constant and 

00 

tj(t) = exp(i7rr/12) TT (1 — exp{i27mr}) (5.65) 

n=l 

is the Dedekind eta function. By combining Eq.s (5.60)-(5.64) and taking into 
account that for e — > + 

r(e) = - £ (l- 1 e + 0(e 2 )) 

and 

tt £ = l + eln7r + 0(e 2 ) 
after some straightforward algebra we obtain (for s — > + ), 

z g( s ) = _1 - sm 2/ [??( T M T )] 2 ) (5.66) 
in agreement with Weil [58], page 75. From here we get 

f s Z- g (s) | s=0 =-ln ? y[77(r)r / (r)] 2 (5.67) 
in agreement with OPS [26]. Eq.(5.9), when combined with Eq.(5.67), produces 

det 'A- g = y\n(T)f. (5.68) 
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Comparison of this result with that given by Eq.s (10.29), (10.30) of Rcf.[20] 
indicates that the free energy T for the Gaussian model on the torus is given 
by 

T =Uny\ V {r)\ 2 . (5.69) 

For the hexagonal lattice A OSP found numerically that det ' A g = ^ \r](j)\ 4 
=0.3 5 5 7 5 23 . In Appendix C it is shown that the hexagonal lattice provides the 
highest packing fraction: q^" = 0.9069. We anticipate that there is a correlation 
between the packing fraction and the numerical value of the determinant (or, 
alternatively, the value of the free energy). For this to happen, according to 
Eq.(5.40) we must have 

det 'A g < dct 'A s = 0.35575 (5.70) 

for any lattice other than the hexagonal. In terms of the free energy defined by 
Eq.(5.69) this means that the free energy of the hexagonal lattice is the highest 
possible 24 . This makes sense thermodynamically. Indeed, according to the G-L 
theory [6] , at criticality, the symmetry of the system is higher than that in the 
low temperature phase. The free energy should be lower in the low temper- 
ature phase to make such a transition to lower symmetry thermodynamically 
favorable. 25 In the next section we shall obtain analogous 3d results for the 
hexagonal close packed (hep) and the face centered cubic lattices (fee) whose 
symmetry is the highest. Such symmetry is typical for the high temperature 
phase (e.g. see Fig.l) and in accord with group-theoretic classification of suc- 
cessive second order phase transitions as predicted (in part) by Indenbom [59], 
e.g. see Fig. 4 in his work. Further details are provided in Section 6 below. 

The above inequality means that we have to prove that it holds for the 
square, hexagonal and triangular lattices. To accomplish this task, we follow 
the classical paper by Chowla and Selberg [29] (C-S) discussed also in our 
earlier work [50] in connection with the Veneziano amplitudes. To facilitate 
the reader's understanding we note that from the point of view of algebraic 
geometry every torus T 2 can be associated with the projective version of the 
elliptic curve whose standard form is given by [51,55,60] 

y 2 =A{x-e{){x-e2){x-e 3 ). (5.71) 

Its periods ui\ and lo-i are given by 

23 We shall recalculate this result below using different methods. 

24 After Eq.(5.3) we have noticed that the free energy is defined with respect to some refer- 
ence state. To be in accord with the G-L theory we can choose as the reference state the state 
at criticality. In this case the value of the free energy for the hexagonal lattice should be taken 
as zero so that the square lattice will have lower free energy typical for the lower temperature 
phase in accord with the G-L theory. Inequality (5.70) indicates that transformation from 
the higher symmetry hexagonal lattice to the lower symmetry square lattice is conformal in 
accord with OSP [26]. 

25 It is tempting to conjecture that the values of respective determinants are proportional 
to the packing fractions. Unfortunately, this is not the case as we shall demonstrate below. 
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lu 1 = fdxUx- ei)(x- e 2 )(x + ei + e 2 )l ^ = 2 -^( fc = ) =; (5.72a) 



e.3 

^= = I dx[( ei - x)(x - e 2 )(x + ei + e 2 )}-i = 2 ^ lK '^ ; (5.72b) 

V/-1 J V e l - e 3 



where k 2 = ^~^ 3 . For the sake of comparison with the results of C-S it is useful 
to keep in mind that the values of integrals K{k) and K'(k) can be rewritten 
in the alternative (Legendre) form, e.g. 

tt/2 

K ^ = I 1 2 (5-73a) 

J 1 — k A sin iyO 


and 

tt/2 

K\k) = K(k') = [ Jf, 2 , (5.73b) 
J 1 — A; 2 sin ip 


where < k < 1 and k 2 + k' 2 = 1. Nevertheless, in the actual calculations shown 
below it is sufficient to use Eq.s(5.72a), (5.72b) as defining equations for K{k) 
and K'(k). Using these definitions we obtain as well, 

r-^if-S. (5.74) 

Naturally, we are interested only in r belonging to the imaginary quadratic field. 
In our case it is cither i or j. Hence, if r is assigned, it is sufficient to know 
only one period in order to determine another. By means of straightforward 
manipulations with elliptic functions C-S demonstrate that for A(r) = [?7(t)] 24 
the following equality holds 



2K 



7T 



12 



A(r)=( — ) 2- 8 (kk'f. (5.75a) 



From here, in view of Eq.(5.70), we obtain 

(A(t))» = (^pj 2-1 {kk'Y . (5.75b) 

In the case of the square lattice the associated elliptic curve is y 2 = x 3 — x 

00 

[60]. This allows us to obtain the period u)\ = J -^0==. Use of substitutions 
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x = 1/y and 1 — y 2 = z, produces w\ = \ J dzz 2(1 — 2;) 4 = I Hi/^llll/ 2 ! . The 

remaining calculations are straightforward and are based on definitions made 
in Eq.s(5.72) and next to it. This allows us to obtain k 2 = k' 2 = 1/2 and, 
accordingly, K (1/^/2) = -j^wi = [r(l/4)] . Collecting terms and using 
Eq.(5.75b) we obtain, 

_L_L T(1 /4V 4 

16 7T 3 



(A(r))^ = --j[r(l/4)] 4 (5.76) 



and, since for the square lattice y = 1, in view of Eq.s(5.70) and (5.75), we 
obtain, det ' A g = 0.3464 . Hence, for the square lattice, inequality Eq.(5.70) 
indeed holds. The OPS result, Eq.(5.70), is given for the hexagonal lattice 
whose dual lattice A* is triangular. Hence, the result, Eq.(5.70), represents 
the determinant for the hexagonal lattice, while the actual calculations were 
made using the reciprocal (dual) lattice, which is triangular. For the triangular 
lattice the associated elliptic curve is y 2 = x 3 — 1 [60]. Therefore, one finds 

k 2 = 5(1 — V — 3) and, because k' 2 = 1 — fc 2 , one finds kk' = 1, which is used 

00 

in Eq.(5.75). Also, for the period uj\ one obtains io\ — J -J===. Substitutions 



x = 1/y and 1 — y 3 = z produce uj\ = ^ / dzz 2 1 (1 — z)e . Accordingly. 



Using Eq. (5.75b) and, taking into account that for the triangular lattice y = 
we obtain 

\3 1 1 \3 rr " 1 ' r ' 




det 'A 3 



2 36tt 2 2 4 / 3 



r(i/6)r(i/2) 



r(2/3) 



0.357567, (5.78) 



to be compared with Eq.(5.70) 26 

The packing fraction for the square lattice (e.g., see Appendix c) can 
be obtained using simple arguments presented in Rcf.[61] and is given by: 

= 0.7853. Accordingly, the ratio /q^ = 1.154845 while the ratio of de- 
terminants is = 1.032237. From here we conclude that, although the 
inequalities between the packing fractions and the inequalities between the re- 
spective determinants are in accord with each other, one cannot claim that the 
values of determinants are proportional to their respective packing fractions. 
Also, there is no need to recalculate determinants for the hexagonal lattice be- 
cause the r parameter is the same as for the triangular lattice [53] and so is 
its imaginary part. We can, loosely speaking, call such lattices self-dual 27 . We 



26 Apparently, there is a minor typographical error in the bound, Eq.(5.70), taken from the 
OPS paper [26]. 

27 According to the accepted terminology [52], the lattice is self-dual if its disciminant is 
equal to l.The cubic lattice in Z" is trivially self-dual (type I lattice). But there are other 
less trivial examples in higher dimensions (type II lattices in terminology of Rcf.[52] ). 
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shall encounter the same kind of self-duality in 3 dimensions as well. This will 
be discussed in the next section. 

6 Ginzburg-Landau theories in dimensions 3 and 
higher 

6.1 Designing higher dimensional CFT(s) 
6.1.1 General Remarks 

In the previous section wc explained a delicate iterrelationship between the 
path integrals Eq.(5.5) and (5.23). From the literature on CFT cited earlier it 
is known that, actually, both are used for desingning of different CFT models 
in 2 dimensions. For instance, if one entirely ignores the effects of curvature in 
Eq.(5.5), then one ends up with the Gaussian-type path integral whose calcu- 
lation for the flat torus is discussed in detail in Ref.[20], pages 340-343, and, 
by different methods, in our Section 5. By making appropriate changes to the 
boundary conditions (or, equivalently, by considering appropriately chosen lin- 
ear combinations of modular invariants) it is possible to build partition functions 
for all exiting CFT models. For the same purpose one can use the path integral 
given by Eq.(5.23), but the calculation proceeds differently as we explained in 
Section 5. Since in 2 dimensions conformal invariance is crucial in obtaining 
exact results, use of the Gaussian-type path integrals is, strictly speaking, not 
permissible. Fortunately, saddle point-type calculations made for the path inte- 
gral, Eq.(5.23), produce the same results since the extremal metrics happens to 
be flat. If one does not neglect curvature effects in Eq.(5.5), one ends up with 
the integrand of the path integral, Eq.(5.23). This result is a consequence of 
Eq.(5.6) known as a conformal anomaly. If one would like to proceed in analo- 
gous fashion in dimensions higher than two one finds that there is a profound 
difference between calculations done in odd and even dimensions. We would like 
to explain this circumstance in some detail now. By doing so we shall provide 
a positive solution to the 2nd Problem formulated in the Introduction. 

In 2 dimensions the conformal invariance of the action, Eq.(5.4), has been 
assured by the transformational properties of the 2 dimensional Laplacian given 
by Eq.(5.1). In higher dimensions, the Laplacian is transformed according to 
Eq.(2.4), so that even the simplest Gaussian model is not conformally invariant! 
This observation makes use of traditional string-theoretic methods in higher 
dimensions problematic. In two dimensions these are based on a two stage 
process. First one calculates the path integral, Eq.(5.5), exactly and, second, 
one uses the result of that calculation (the conformal anomaly) as an input 
in another path integral, e.g. Eq.(5.23), which is obtained by integrating this 
input over all members of the conformal class. Since in odd dimensions there is 
no conformal anomaly as we shall demonstrate momentarily, such a two stage 
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process cannot be used. Moreover, absence of a conformal anomaly in odd 
dimensions also affects the results of finite scaling analysis 28 . The situation 
can be improved considerably if we do not rely on the two stage process just 
described. We would like to explain this fact in some detail now. 

Even though the transformational properties (with respect to conformal 
transformations) of the Laplacian, Eq.(2.4), in dimensions higher than 2 are 
rather unpleasant, they can be considerably improved if, instead of the usual 
Laplacian, one uses the conformal (Yamabe) Laplacian O g defined by 

a g = A g + aR(g) (6.1) 

where a = aT 1 = 4^— ?• By construction, in 2 dimensions it becomes the 

4 d—1 J 1 

usual Laplacian. In higher dimensions its transformational properties are much 
simpler than those for the usual Laplacian (e.g see Eq.(2.4)). Indeed, it can be 
shown [62] that 

D e 2 fg = e-^+i)/n se (l-i)/. (6.2a) 

This result can be easily understood if we use results of Section 2. Indeed, since 
e 2 ^ = tp p ~ 2 and p = ^4j, we obtain at once e^^ 1 )/ = ^, while the factor 

e -(i + 1 ) is transformed into Lp x ~ p . From here, it is clear that Eq.(2.12) for scalar 
curvature can be cquivalently rewritten as 

R(g) = a^-^Agcp + a- l R(g)ip), (6.3) 

where g = e 2 ^ g , so that 

R(g) - aU e 2i g . (6.4) 

In practical applications it could be more useful to consider two successive con- 
formal transformations made with conformal factors e 2 * and e 2h . If e 2 ^ — Lp p ~ 2 
and e 2h — i}j p ~ 2 then, we obtain, 

A§V + or x R{g)il> = ^ p [A 3 (W) + oT x R{g) (W)]- (6.2b) 

This result is, of course, equivalent to Eq.(6.2a). 
Consider now the following path integral 

exp (-H9)) = J D[p] C MSu g &)} (6.5) 

where 

S D » = f d?x y /g{(v g <p) 2 +cr 1 R{<,)<p 1 } 
Jm 

= f d d x^g[ipn g cp] = f d d x^g~R(g). (6.6) 
Jm Jm 

The conformal factor (p~ p in Eq.(6.3) is eliminated by the corresponding factor 
coming from the volume factor of y/g = ^/exp(2/)g = (p P y/g, as explained after 

28 As mentioned in Section 5.2, even in 2 dimensions the results of calcultations of conformal 
anomalies were used incorrectly in the finite size scaling analysis. 
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Eq.(3.5). Thus, Eq.(6.5) is the exact higher dimensional analog of the two di- 
mensional path integral, Eq.(5.5). Thus, problems related to higher dimensional 
CFT are those of Riemannian (quantum) gravity [25]. In this paper we are not 
considering the pseudo Riemannian case associated with Einsteinian gravity. 

The question arises: If the path integral, Eq.(6.5), is such an analog, is there 
a higher dimensional analog of Eq.(5.6)? The answer is "yes", if the dimension 
of space is even and " no" if the dimension of space is odd [62] . 

Because Eq.(5.6) is used heavily in finite size calculations [44], the absence 
of similar results in 3 dimensions should be taken into account. Earlier attempts 
to generalize these two dimensional results to higher dimensions [44] were made 
without such consideration. 



6.1.2 Lack of conformal anomaly in odd dimensions 



In view of its importance, we would like to provide a sketch of the arguments 
leading to the answer "no" in dimension 3, important for developments in 
this paper. Clearly, the same kind of arguments will be of use in other odd 
dimensions. In doing so, although we follow arguments of Rcfs. [62,63], some of 
our derivations are original. We begin by assuming that there is a one parameter 
family of metrics: g{x) =exp(2x/)g. Next, we define the operator 8 f via 

5fUg = ~dx~ ' x=0 Dox p( 2:e /)9 ' ( 6 ' 7 ) 
Taking into account Eq.(6.2) we obtain explicitly 

6 f D g = -2fU g (6.8) 

and, accordingly, 

5 f e-^ = -t(5 f U g )e-^. (6.9) 
These results allow us to write for the zeta function (Appendix B) 

oo 

*/Cn» = ^ J dU^SfTrie-^) 



o 



'JO 







dtt s Tr(-2fa g e 



-ta„ 



-2s 
o 



J dtt s - l Tr{fe- ta °). (6.10) 



The last line was obtained by performing integration by parts. Since for small 
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t's it is known that, provided that dimkerD g = 29 , 

Tr{fe-^) ~ f(x)u k (x)dv<*J t k ^ 2 . (6.11) 

Using this result in Eq.(6.10) produces 

1 CO 

SfCn g (0) = --^{J dtt^Tr{fe- tD «)+ J dtt^TrUe-^)) 

1 

OO 

- -rW<|fe^V d<i '" ,rr(/e "' D ' ))l " 0<6 ' 12) 

Since for s — ► + we have (l/r(s)) ~ s, the second (regular) term in brackets 
will become zero when multiplied by the combination j^fy, while the first term 
will become zero even if it might acquire a pole (when d = 2k). Hence, for all 
dimensions d > 3 we obtain <5/Cn g (0) = or 

Cn,(0) = C Ds (0). (6-13) 
This result can be used further now. Indeed, if we write 

Sf [rWCn g W] = T(s)[5 f C as (0) + sS f C' a3 (0) + O(s 2 )} (6.14) 

and take into account Eq.(6.13) and the fact that sT(s) = 1 we obtain, 

CO 

<5/C'n s (0) = S f J dtt^Trie-^) | s=0 





oo 



= - 2s (E i + k-d/2 + J dttS ~ lT ^- tDs )) l-o • (6-15) 

Applying the same arguments to Eq.(6.15) as those which were used for Eq.(6.12) 
we conclude that, provided that dimkerD g = 0, in odd dimensions, <5/Cn (0) = 
0. That is 

Cn 9 (0) = Cn 9 (0)- (6.16) 
In view of Eq.(5.9) this leads also to 

detn s = detn g , (6.17) 

QED. 



29 For path integrals this is always assumed since zero modes of the corresponding operators 
are associated with some kind of translational, rotaional, etc. symmetry. To eliminate the 
undesirable dilatational symmetry, one actually should use the Yamabe functional, Eq.(3.6), 
as explained in Section 3. Although this is silently assumed thus far, arguments additional to 
those in Section 3 will be introduced further below. 
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6.1.3 3d CFT path integrals 



The previously obtained results can be refined further if we recall Eq.s(5.9)- 
(5.11). In particular, let e 2 ^ in Eq.(5.11) be rewritten as some nonnegative 
constant I. Then we obtain 

(~ g (s) = l s C g (s). (6.18) 

This result is consistent with Eq.(6.13) for s = 0. Differentiation with respect 
to s produces 

C~(0HC s (0)lnZ + C;(0). (6.19) 

In view of Eq.(5.9), this result is equivalent to lndct D g = lndet D g — 
C g (0) In/. This result apparently contradicts Eq.(6.17), but the contradiction is 
only apparent in view of the earlier footnote. The situation is easily correctable 
if in the path integral, Eq.(6.5), we replace the action functional Sa(ip) by that 
of Yamabe given by Eq.s(3.2) (or (3.6)). This, by the way, allows us to fix the 
value of C 9 (0) m Eq.(6.19): provided that we identify the constant I with the 
volume V, the value of ( g (0) = |. After this, the situation in the present case 
becomes similar to that encountered in the previous section, e.g. see Eq.s (5.12 
and (5.13). Now, instead of the functional T{<£>) given by Eq.(5.15), we consider 
the related functional given by 

T{tp) = lndctDg -C s (0)lnV 

= IndetCL - -lnV. (6.20) 
V 

For the path integral calculations, a functional defined in such a way is not 
yet sufficient. To repair this deficiency we have to impose a volume constraint. 
That is we need to consider the path integral of the type 

Z Y (V) = f D[<p]6{ f d d x^ - V)exp(-S[<p}) (6.21) 

with S[ip] given by Eq.(3.2) (or (3.6)). 

Clearly, the path integral Zy{V) (Y in honor of Yamabe) is the exact higher 
dimensional analog of the "stringy" path integral Sx(-A) given by Eq.(5.33). In 
view of Eq.(3.7), it also can be viewed as the path integral for pure gravity in 
the presence of the cosmological constant. Because of this, the standard path 
integral for the self interacting scalar ip 4 (or LGW) field theory is obtainable 
now in complete analogy with Eq.(5.32), i.e.: 



/ 



OO 

D[4>]e^{~S L Gwm = J dVe- bv Z Y (V) (6.22) 

o 



But, since the variation of the Yamabe functional produces the same Eq.(3.3) 
as can be obtained with LGW functional, Slgw{<I>), one can develop things 
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differently, but surely equivalently. To this purpose, instead of the functional 
S[ip] given in Eq.(3.2) we use 

SvM = X [ d d x^g{( Vglf f + R(g)cp 2 } (6.23) 

and replace S[ip] in the exponent of the path integral in Eq.(6.21) by Sv [f] from 
Eq.(6.23). Then, instead of Eq. (6.22), we obtain, 

OO 

J D[4,]e X p{-S LGW (<P)}= J dVZ Y (V), (6.24) 

o 

where the sign = means "supposedly" . This is so, because, at the level of saddle 
point calculations the left hand side and the right hand side produce the same G- 
L equation. Beyond the saddle point, calculations are not necessarily the same. 
Although we plan to discuss this issue in detail in subsequent publications, some 
special cases are further discussed below in this section. 

It should be clear, that at the level of saddle point calculations, replace- 
ment of the functional S[ip] in Eq.(6.21) by from Eq.(6.20) is completely 
adequate, so that the sequence of steps in analysis performed for the two dimen- 
sional case in Section 5 are transferable to higher dimensions without change. 
This can be summarized as follows. 

Although in 3 dimensions we have the result given by Eq.(6.17), which 
forbids use of identities like that in Eq.(5.6), still, based on arguments just 
presented, the functional T((p) defined by Eq.(6.20) should be used in the ex- 
ponent of the corresponding path integral replacing that given in Eq.(5.15) in 2 
dimensions. Since by doing so one will be confronted with the same type of min- 
imization problems as discussed earlier in Section 5 30 , thus defined functional 
integral is an exact 3 dimensional analog of the path integral, Eq.(5.30). 



6.2 Completion of the work by Lifshitz in 3 dimensions 
6.2.1 General remarks 

The results just obtained allow us to proceed with the rest of our developments 
in complete accord with results of Section 5. This means that our task from 
now on will be to provide an affirmative answer to the 1st Problem formulated 
in Section 1.4 following logical steps developed previously. To this purpose, 
it is helpful to use some results from the classical paper by Hawking on zeta 
function regularization of path integrals in curved spacetimc [65] in view of the 
fact that, typically, the calculation of path integrals is done by the saddle point 

30 This conclusion had been reached without any reliance on path integrals and on physical 
applications in Rcf. [64]. In this and related Rcf. [31] the extremal properties of determinants 
of O g and A g with respect to variations of the background metric were studied. 
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method. In the present case, the saddle point level of approximation for the 
path integral, Eq.(6.21), is equivalent to minimization of the functional f(ip) 
given by Eq.(6.20). In turn, this is equivalent to minimization of the Yamabe 
functional, Eq.(3.2). The following theorem (attributed to Aubin, Trudingcr 
and Yamabe) is very helpful for this task. 

For any d- dimensional compact manifold M 

A(M) < X(S d ) (6.25) 

where X(S d ) is the Yamabe invariant for d-dimensional sphere. 
The proof of this result can be found in Refs. [22,23]. 

The previously obtained inequality, Eq.(5.37), indicates that, at least in 
two dimensions, the same type of inequality exists for determinants. These 
two dimensional results for determinants cannot be readily extended to higher 
dimensions however since, in view of Eq.(6.17), there is no conformal anomaly 
in 3 dimensions (but there is in dimension 4 [66,67]!). Therefore, in Refs. 
[66,67] extremal properties of these determinants with respect to changes in 
the background metric have been studied. These studies produced inequalities 
analogous to that given in Eq.(6.25), which take place only locally, i.e. in the 
vicinity of some point onp€ M . 

The Yamabe constant X{S d ) has been calculated. In particular, for S 3 it is 
found to be X(S 3 ) = 6(27r 2 )i, Rcf.[33] 31 . Calculation of X(M) for manifolds 
other than those diffcomorphic to S 3 has been found to be surprisingly difficult. 
Only few explicit examples are known to date [33]. Fortunately for us, they are 
just sufficient for our current purposes. In particular, Gromov and Lawson [68] 
have demonstrated that X(T d ) = for d > 3 and, moreover, the same result 
holds for the connected sum (T d # T d , etc.) of d— dimensional tori. Aubin 
Ref.[22], pages 150-152, proved the following theorem 

Let M be d-dimensional C°° compact Riemannian manifold, then there is a 
conformal metric whose scalar curvature is either a nonpositive constant or is 
everywhere positive. 

From the results of Section 1 and 2 it should be clear that in the case of the 
G-L functional, based on physical arguments (e.g. read the discussion related to 
Eq.(1.20)), one should look at cases of nonpositive constant scalar curvatures. 
That is, one should look at the mass term m 2 representing constant negative 
scalar curvature for temperatures below criticality (and zero scalar curvature at 
criticality) in accord with Eq.(2.10). 32 . 

Using the above Theorem by Aubin and taking into account details of its 
proof we conclude that for manifolds admitting constant curvature (spherical(s), 
hyperbolic(Zi) or flat (0)) one should expect 

2 

31 In d-dimensions this result is replaced by \(M) < d(d — 1) [yoi(5 d )] d . 

32 In our recent work, Ref. [69], we have studied the connections between the hyperbolicity 
and conformal invariance. Development of such connections is linked with the study of prop- 
erties of hyperbolic 3-manifolds and orbifolds. Examples of such connections can be found in 
the same reference in the context of AdS-CFT correspondence. 
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(6.26) 



This result is consistent with that known in two dimensions. Indeed, formally 
using Eq.(3.6) in two dimensions and employing the Gauss-Bonnet theorem we 
obtain 

A(M) = Attx(M). (6.27) 

As discussed in Section 1 in this work, we are interested mainly in deter- 
mining symmetries of the high temperature phase. In this phase the G-L order 
parameter is zero and so is the free energy at the saddle point level of approxi- 
mation (e.g. read the comments after Eq.(5.31)). This implies that the Yamabe 
constant is zero, i.e. that X(T d ) = 0. Next, using Eq.(6.16), valid only if the 
volume constraint is imposed as we explained above, for a given background 
metric g we can choose the metric g as flat. With such a choice of metric we 
can use the results of Richardson, Ref.[31], which we would like to summarize 
briefly now. 

Let g(u) = e^ u ^g be 1-parameter family of metrics of fixed volume and such 
that .9(0) = g. This implies that 0(0) = and / e^dVo = V. In two dimcn- 

M 

sions, using results of OPS [26], especially their Eq.(1.12), it is straightforward 
to obtain the following result 

A(-lndct > Ag(u) ) | u=0 = C>) = T^/ j> K (9)dV g , (6.28) 

M 

where K (g) is the Gaussian curvature for the metric g. C g (0) and <j> represent 
SlCg(t I )(0) U=o an d ^0( u ) | u =orespectively. Volume conservation implies 

|- Je^dV g | u=0 = f Wg = (6.29) 

M M 

in accord with the earlier result, Eq.(5.17). If the Gaussian curvature K(g) is 
constant, then Eq.(6.28) and Eq.(6.29) produce the same result. This implies 
that J^C s («)(0) |u=o= 0. That is g, is the "critical" (extremal) metric. In view 
of Eq.(6.28), this also means that for such a metric the free energy attains its ex- 
tremum. We encountered such extremal situations in the previous section when 
we demonstrated that the triangular and hexagonal lattices have the highest 
possible free energies as compared to other lattices. We shall demonstrate be- 
low that the same holds in 3 dimensions where the hep and fee lattices possess 
the highest possible free energies as compared to other lattices. To this purpose, 
we need to generalize the above two dimensional results to higher dimensions. 
Formally, this is not an easy task mainly because of the differences in properties 
of the Laplacian under conformal transformations in two and higher dimensions, 
e.g. see Eq.(2.4), causing absence of the Moser-Trudinger type of inequalities in 
dimensions higher than two. Nevertheless, some of the difficulties can be easily 
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overcome with help of the results we have already. For instance, by combining 
Eq.s (6.10), (6.14) and (6.15) we obtain at once 

oo 

<5/Cn 9 (0) = S f J dtt^Trie-^) \ s=0 



oo 

= J dtt s Tr(-2fn g e- t0 °) | s=0 . (6.30) 
o 

From here we obtain essentially the same result as the main theorem by Richard- 
son [31], i.e. 

C' Ds (0) |„= = = J dV g j>(x)n g ((l, x, x). (6.31) 

M 

That is, provided that we replace \3 g by A g and require that locally £(1, x, x) = 
const, with the heat kernel ((s,x,x) given (as usual) by 

C( W ) = £4P' ( 6 - 32 ) 
fc=i Afe 

we obtain the main result by Richardson, Ref.[31], e.g. see his Theorem 1 
and Corollary 1.1. Here ip k (x) are eigenfunctions of the Laplacian (or Yamabe 
Laplacian respectively) corresponding to eigenvalues A^ with = A < Ai < 
A 2 < • • • < Afc < • • •. In two dimensional case the condition for criticality is given 
by Cg(0) = is local meaning that, provided the the volume is constrained, the 
constancy of the Gaussian curvature K(g) at given point of M is caused by the 
metric g for which C ff (0) is extremal. In 3 and higher dimensions, the constancy 
of curvature at the point of M is replaced by constancy of £(l,x,x) under the 
same conditions of volume conservation. This condition is necessary but is not 
sufficient now, since the analog of the Moser-Trudinger inequality (used to prove 
sufficiency in 2 dimensions) does not exists. Instead, one should study locally the 
second variation of C' g (0) with respect to the underlying background metric in 
order to decide if such (local) extremum is maximum or minimum. Fortunately, 
this task was accomplished in Ref.s [31,64]. In particular, Richardson [31] 
obtained the following theorem of major importance for our work: 

The Euclidean metric on a cubic 3-torus is a local maximum of determinant 
of the Laplacian with respect to fixed- volume conformal variations of the metric. 

This Theorem is proven only for the cubic 3-torus. The word "local" means 
that there could be (or, there are, as we shall demonstrate) other 3-tori also 
providing local maxima for determinants. In fact, according to the result by 
Chiu [70], all determinants of flat 3-tori possess local maxima so that the de- 
terminant for the face centered cubic (fee) lattice has the largest determinant. 
Unfortunately, his results are nonconstructive and hence, cannot be used in 
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physical applications. Therefore, below we provide an entirely different way to 
reach the same conclusions. 

The second variation of the Yamabe functional was calculated by Muto [71] 
(see also Ref.[25]) with the result: 

M 

where the constant a = d — 1. As in the case of quadratic actions in the flat 
space [15,19] the second variation (with volume constrained to be equal to one) 
looks very much the same as the original quadratic Yamabe functional, except 
for the "wrong" sign in front of scalar curvature 33 . Following Muto [71] we 
conclude that: a) if R{g) is positive, the second variation can be made positive 
for appropriately chosen <p 34 , b) if R(g) is negative, it is positive for the same 
reasons. 

The positivity of second variation implies that the extremal constant cur- 
vature metric g provides a locally stable minimum for lZ(g(t)) (that is, using 
results of Section 4, the Einstein metric obtained as solution to Eq.(4.2) is 
stable among nearby metrics). 

It is interesting to notice that calculation of higher order fluctuation correc- 
tions to the Yamabe path integral, Eq.(6.21), involves calculations on the moduli 
space of Einstcinian metrics, Ref. [72]. This observation provides a strong link 
between higher dimensional LGW theory and two dimensional string inspired 
CFT discussed in the previous section. Naturally, Eq.(6.24) can be used to inves- 
tigate to what extent the final results of conventional field-theoretic calculations, 
Ref. [19], may differ from more sophisticated string-theoretic calculations in the 
style of Ref.s [27,28,72]. This task is left for further study. 

For the same reasons as in the two dimensional case considered in the previ- 
ous section, we are not interested in a positive constant curvature metric typical 
for Steven though such metric do provide a local maximum for determinants 
of both Laplace and Yamabe operators [64]. Thus, we are left only with zero 
and negative curvature metrics characteristic for physical systems at and below 
criticality (i.e. below T c ). Thermodynamically, the system at criticality should 
possess higher free energy than the system below criticality. In our case, at crit- 
icality the free energy T=\ lndet (e.g. see Appendix D). Hence, the above 
cited theorem by Richardson guarantees that for cubic 3-tori the free energy 
is larger as compared to that below criticality (which corresponds to negative 

33 This "wrong" sign has significance, however. It is in accord with existing calculations of 
the fluctuation corrections to G-L theory in the low temperature phase, e.g. sec Eqs.(1.20), 
(1.22). Moreover, given that the low temperature "mass" term in Eq.(1.22) is 2 \a\ and 
comparing the expansion, Eq.(1.21), with Eq.(6.33) we obtain (for d = 4) exactly the same 
fluctuation kernel using Eq.(6.33). This is so since in d = 4 the combination ^^-R(g) = 2R(g) 
and R(g) is negative in the low temperature phase. The constant c in front of Laplacian is 6 
in the present case while it was left unspecified in Eq.(1.22). 

34 This can be easily understood if we expand ip into Fourier series made of eigenfunctions of 
the Laplacian and take into account that for any closed Ricmannian manifold the spectrum 
of the Laplacian is nonncgative and nondecreasing [22]. 
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curvatures) and Chiu's results [70] imply (by analogy with two dimensions) 
that among 3-torus lattices the fee lattice possess the highest free energy. Be- 
low we shall provide more direct demonstration of this fact and, in addition, 
we shall demonstrate that the spectra of both fee and hep lattices possess the 
same value for determinants, thus implying the possibility of phase transitions 
between these lattices. Such transitions have been indeed observed in nature 
[73], but their description falls outside of the scope of this paper. 

Thus, in the rest of this section we shall concentrate on explicit calculations 
of the toral determinants for different 3-dimensional lattices. Such calculations 
involve the use of the 3-dimcnsional version of the zeta function, Eq.(5.41). 
Unlike the two dimensional case considered earlier, we cannot apply directly 
the 1st Kronecker limit formula in order to obtain the corresponding values for 
determinants. We are also unable to extend the ideas of complex multiplication 
in order to determine the types of allowable lattices. For even dimensional 
spaces, where one can use the concept of the Abelian variety such a task can be 
accomplished. Some of these varieties possess complex multiplication and can 
be mapped into even dimensional tori [74] . We are not aware of similar results 
for odd dimensional spaces. Thus, we have calculated determinants numerically 
using a procedure to to be discussed below. The results obtained are in complete 
qualitative agreement with those obtained in two dimensions. 

In order to introduce our readers to issues involved in such calculation, we 
would like to discuss the 3 dimensional analog of the 1st Kronecker limit formula 
now. It was considered by Bump and Goldficld [57] and later summarized in 
the paper by Efrat [75] . 

We begin with the following observation. As it is shown by Sarnak [76] , the 
function E(t,s) in Eq.(5.61) is an eigenfunetion of the hyperbolic two dimen- 

d 2 d 2 

sional Laplacian A = — w 2 (-r— ~ + tt-^), i.e. 

ox z oy z 

AE(t,s) = s(l- s)E(t,s). (6.34) 

We discussed this type of eigenvalue equation extensively in our earlier work, 
Ref.[69]. In the same work we discussed the d— dimensional extension of such an 
eigenvalue problem. Unfortunately, these results cannot be used in the present 
case. As it is demonstrated by Bump [77], instead of looking for eigenfunc- 
tions of the 3 dimensional hyperbolic Laplacian, one should consider a more 
complicated eigenvalue problem. It is important to realize at this point that 
this eigenvalue problem should be of the same relevance to all exactly solvable 
3 dimensional statistical mechanics models as two dimensional eigenvalue prob- 
lem, Eq.(6.34), to two dimensional exactly solvable models discussed in Rcf. [20]. 
Since we are not aware of exact solutions of non trivial 3 dimensional models, we 
believe that discussing the issues involved in such calculations might shed some 
new light on the whole problem of exact solvability in dimension 3. In addition, 
we would like to present these results in order to compare them against our 
calculations of determinants presented below. 

To begin, we need to write the 3 dimensional analog of E(t, s) (e.g. compare 
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with Eq.(5.61)). It is given by [78] 



E(r,s) = ]T — — M . (6.35) 

(m,n,fe)=i \ kz2 + m \ + i kx 3 + mxi + n) 2 } 2 

In Appendix B we introduce the Epstein zeta function, Eq.(B.5). It can be 
demonstrated [57,75] that 

Z(A,s) = |Af f C(3s)E(t,s) (6.36) 

where |A| = det A and Z2 = X2 + iyi, Hi, 1/2 € R>o, x\, X2, Xs € R. By anal- 
ogy with the two dimensional case, one can establish a very important func- 
tional equation for the combination E*{t 1 s) = tt^~T(^-)((3s)E(t, s). Namely, 
E*(t,s) = E*(t, 1 — s) [57]. This equation allows us to extract the s — > 0+ 
limit for E(t, s) using the associated result for s — > 1 + obtained in Efrat's paper 
[75]. For the combination ((3s)E(t, s), which he also denotes as E*(t,s), he 
obtains 

= + (C - 2/3 In ( K yl) - 4 taj(r)) + 0(« - 1), (6.37) 



where 

1/2 

g(r) = exp(— 1 2 E*(zi,s)) H |1 — cx p(— 2iry2 \nz\ — m\ + 2ni(mxi + nx^})\ 

OTT 

(6.38) 

with z\ = X\ + iyi and X4 = X3 — X1X2 and C being a known constant. The 
limiting expression (for s — > + ) can be found in the work by Chiu [70]. Since 
neither Chiu nor Efrat have provided any explicit examples of actual calculations 
involving formulas just presented, we have chosen another approach to the whole 
calculation of these limits. Before discussing our calculations we would like to 
mention that Bump [77] and, following him, Efrat [75] have demonstrated that 
E(t,s) defined in Eq.(6.35) is an eigenfunction of two operators Ai and A2 
(whose explicit form is rather complicated [77]) so that 

AiE(t,s) = 3s(s-1)E(t,s) (6.39a) 
A 2 £(t,s) = -s(s-1)(2s-1)E(t,s) (6.39b) 

to be compared with Eq.(6.34). 



6.2.2 Calculation of determinants: general discussion of the numer- 
ical algorithm 

So far in this paper calculation of determinants has been done with the 
help of zeta functions, e.g. see Eq.(5.9). This method required use of the first 
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Kronecker limit formula in two dimensional calculations. Generalization of this 
result to higher dimensions was formally accomplished by Epstein at the turn of 
20th century [79]. However, his results were too general to allow any practical 
calculations. This fact caused many attempts to improve/simplify his calcu- 
lations. The result, Eq.(6.37), is one of many [70,80], etc. It can be derived 
directly from the 3 dimensional Epstein zeta function [75,77]. However, sub- 
sequent results by Chiu [70] stop short of actual use of Eq.(6.37) in order to 
produce numerical values for lattice determinants, as it is done in two dimen- 
sions. In view of this, we are going to develop another method of calculation of 
lattice determinants, which does not involve the use of zeta function. The relia- 
bility of this alternative method is tested against exactly known results obtained, 
again, with the help of some ideas from physics. In Appendix D we provide the 
simplest example of this type of calculation, so that readers are encouraged to 
read this appendix before proccding with the rest of this subsection. 

Assuming this, we need to calculate the 3 dimensional version of Eq.(5.41), 

i.e. 

Zg(s)=y^' r-g . (6.40) 

9 ^' £A * (47T2 |1| 2 ) 

As in Appendix D, we disregard the constant An 2 in the denominator and con- 
sider the following regularized sum instead 



V |1*| +K 2 



(6.41) 



where now the summation takes place over all lattice cites of the reciprocal 
lattice A* due to the presence of the small parameter k 2 , which will be put to 
zero at the end of calculations. Evidently, as in Appendix D, we have 



(3T = lim - In 



l,_detA 5 (K 2 ) 



K 



o2 detA § (0) J ^\\*\ 2 + k 
o i 1 1 

Again, as in the Appendix D, to proceed, we need to evaluate somehow the sum 
at the r.h.s. of Eq.(6.42). To this purpose we would like to take advantage of 
the Poisson summation formula [81]. In the present case we obtain, 



V- 



1 / vol{A) ^ e-*! 1 ! 



where we take into account that in 3 dimensions the summand on the l.h.s. rep- 
resents the Fourier transform of the screened Coulomb (also known as the Debye- 
Hiickel (D-H)) potential displayed as the summand on the r.h.s. The parameter 
k is known in the literature on electrolyte solutions [6,82] as the D-H inverse 
screening length. Treated from such perspective, our calculations are reminis- 
cent of those for the Madclung constant in solid state physics [13, 83]. Using 
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Eq.(6.43) in Eq.(6.42) and taking into account results of the Appendix A we 
obtain upon integration over k the following result: 



(& = lim(-l) Kv ^(A) £ &M J 6 ' 10 - ( 6 - 44 ) 

1 W 

This result is still very inconvenient to use. To simplify matters further we make 
use of the expansion Eq.(B.9) for the theta function and take also into account 
Eq.s(B.5)-(B.8). This produces the following result to be used in numerical 
calculations: 

QT= Um(-l)K A /disc(A)^ — exp(-«V^). (6.45) 

n=l 

As in Eq.(B.7), we have subtracted one so that the series starts with n = 1. The 
coefficients a(n) have been tabulated for various lattices [52, 84]. The difficulty 
in evaluating such sums lies in the fact that there is only a finite number of a(n)'s 
which are available in literature. The difficulty with calculations containing 
finite number of terms can be seen already in evaluation of much simpler sums 
such as, for example, 

oo 

S(k) = k cxp(— Kfi). (6.46) 

n=0 

Clearly, evaluating the geometric progression and taking the limit: k — > + 
produces 1 as expected but, should we keep a finite number of terms and take 
the same limit, we would obtain an entirely wrong result. The results can be 
considerably improved if we correlate the number of terms in the sum N with 
the optimal value of k for such N. This can be achieved by minimizing the 
above sum with respect to n. This leads to the following minimization equation: 

JV N 

cxp(-«*n) = k* ^2 n cxp(-K*n) (6.47a) 

71=0 71=0 

or, equivalently, 

\ = (n). (6.47b) 

The best results are obtained with still additional refinement. Since we are 
interested in the limit k — > + , it is convenient to look at the obtained 
values of k* as function of N, i.e. we look for k*(N). We expect that for 
N 2 > Ni the optimal n's should behave as k*(N 2 ) < k*(N\), etc. This indeed 
happens. Moreover, we expect that S(k*(N)) behaves in a similar way, i.e. 
if we have successive sums including N\ < N 2 < N$ < • • •, we expect that 
|S(k*(JVi)) - S(k*(N 2 ))\ > \S(k*(N 2 )) - S(k*(N 3 ))\ > • • • so that the series 
is converging in the desired direction. For the sum like that given in Eq.(6.46) 
this is easy to check, but for the sum, Eq.(6.45), this is less obvious (although 
numerically we obtained exactly the same type of convergence). So, additional 
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arguments need to be invoked. For instance, the sum in Eq.(6.44) can be an- 
alyzed as is typically done in solid state physics. Here, one routinely replaces 
such sums by integrals [13,83]. In our case we obtain, 

oo 

Eexp(— k|1|) f , ,exp(— kx) ,„ 
1|2 ~ « J <*™ 2 X 2 ~ !■ ( 6 - 48 ) 
ill o 

As plausible as it is, one can obtain a much better estimate based on the work 
by Jones and Ingham [85]. It is more accurate when used for the sum, Eq.(6.45). 
According to these authors, if one thinks about a selected atom at the lattice site 
chosen as the origin, then at large distances R from the origin, it is permissible 
to assume that the atoms are more or less homogeneously distributed on the 
sphere of radius R. In view of this, the factor a(n) is proportional to the area 
of a sphere of radius R while the factor of n <~ R 2 . Hence, in our case the sum 

oo 

S{k) =k^2 exp(-KVn) (6.49) 

71=1 

provides an upper bound for the sum in Eq. (6. 45). The lower bound can be 
estimated if we replace all a(n)'s in Eq.(6.45) by 1 thus obtaining the sum 

S(k) — k — exp(— Ky/n). (6.50) 
n 

n=l 

6.2.3 Use of designed algorithm for calculation of various lattice 
determinants 

The convergence procedure for the sums defined above can be easily tested 
for any N and, indeed, it satisfies the convergence criteria just described. A 
numerical check provided assurance that the sum, Eq.(6.45), is indeed bounded 
by the above two sums for lattices we are interested in. The numerical values for 
coefficients a(n)'s are taken from Rcf.[84]. Specifically, Table 7 of this reference 
provides 80 entries for the cubic lattice, Table 11 provides 80 entries for the fee 
lattice, Table 16 provides 80 entries for the hep lattice and Table 24 provides 80 
entries for the bec lattice. The book by Conway and Sloane (C-S), Ref. [52], 
provides the values of discriminants for lattices of unit edge length. In particular, 
for the bec lattice we find: 'disc(A) = 4, and accordingly 1 for the cubic, 2 
for the fee and \[2 for the hep lattices. In actual computations the following 
information was taken into account. First, both the cubic and hep lattices are 
self-dual [52]. Second, the bec lattice is dual to fee [13,52]. This means, for 
instance, that one should be careful when using the relation vol(A)vol(A*) = 1 
(or, as in solid state physics, Ref. [13], vol(A)vol(A*) = (2tt) ) since, if one 
chooses volAbcc = 4 then, one formally gets volA*j cc = 1/4 (or (2ir) /A) as 
compared with 2 for fee just presented. To resolve this "paradox" we use results 
from the solid state physics [13]. For instance, the basis vectors of the fee lattice 
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are known to be: ai = §(— a x + a z ), a 2 = §(a y + a z ),a 3 = §(— a x +a y ). From 

here the volume is obtained as V/ cc =ai • (a 2 x a 3 ) = Using this result for 
a = 1, we obtain: U/ cc = 1/Vb cc where Vb cc = 4 in accord with C-S. At the 
same time, C-S provide the value 2 for Vf cc . As a final example, let us consider 
the bcc lattice. In this case we have the basis vectors ai = | (— a x — a y + a z ), 
a 2 = %(a y + +a y + a z ), a 3 = | (— a x + a y — a z ) so that the volume is obtained 
as Vbcc = if- From here, we obtain V& cc = 1/V/ CC where Vf cc = 2 in accord 
with C-S. Since the hep lattice is selfdual [52], we obtain Vhcp = l/V/j Cp = 
1/V2- These examples illustrate the relationships between the (C-S) and the 
accepted solid state (SS) conventions. They are summarized in the Table 2 
below. 



Table 2 



Lattice 


Vol 


Laticc 


Vol 


bcc{C - S) 


4 


fcc*(SS) 


1/4 


fcc(C - S) 


2 


bcc*(SS) 


1/2 


Z^C-S) 


1 


Z*(SS) 


1 


hcp(C - S) 


V2 


hcp*(SS) 


1/V2 



It should be clear from this table and the examples just presented that 
the conventional solid state direct lattice results are dual to that given in C-S 
monograph [52]. To make sense out of our numerical results, we use the C- 
S results in our calculations compatible with C-S data for a(n)'s. Using the 
numerical procedure outlined above the following results for the free energies 
are summarized in Table 3 below. 



Table 3 



LatticeA* 


Free energy 


Packing fraction 


hep 


-6.669431 


0.74 


fee 


-6.65616 


0.74 


bcc 


-8.68137 


0.68 




-9.27008 


0.52 



These results are in qualitative agreement with those obtained earlier in two 
dimensions where we observed that the respective free energies are arranged in 
accordance with their respective packing fractions. The accuracy of our calcula- 
tions is supported by the fact that we have obtained numerically practically the 
same free energies for the hep and fee lattices whose packing fractions are the 
same (e.g. see Appendix C). As in two dimensions, the earlier made assump- 
tion that the free energies are proportional to their respective packing fractions 
happens to be wrong. The obtained results provide an affirmative answer to the 
1st Problem formulated in Section 1.4. 

The equality of free energies between the hep and fee lattices has been 
established in recent computer simulation methods [86] and apparently is of 
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great practical interest. Such an equality between the free energies might lead 
to the phase transition between these lattices and, indeed, they were observed 
and theoretically discussed [73] . 

The equality between the hep and fee free energies has been achieved with 
account of the following observation. On page 114 of C-S book, Ref.[52] it is said 
that: "The hep is not itself a lattice, but may be defined as a union of a lattice 
L....and the translate...". Such an opinion is not shared by solid state physicists 
for whom hep is a legitimate lattice [8,13,83]. If one compares between C-S 
and SS definitions, then one should multiply \f2 in Table 2 by a factor of 2. 
After such multiplication the free energies of both hep and fee lattices become 
the same. The data of Table 3 reflect this observation and are in accord with 
experimental results depicted in Fig.l, which we would like to discuss now. 

6.2.4 6.3. A brief walk across the CuZn phase diagram. 

As is mentioned in the Introduction, the CuZn phase diagram is a typical equilib- 
rium phase diagram for a binary alloy. Many such alloys make superlattices un- 
der appropriate concentration/temperature conditions[8,87]. Apparently, this 
fact was a motivation for Lifshitz original work, Rcf. [4]. Taking this into ac- 
count, we would like to complete this section by connecting the phase diagram 
depicted in Fig.l with the results summarized in Table 3. 

To do so, we would like to remind our readers about how such phase dia- 
grams are obtained in real life and what they are actually supposed to convey. 
Initially, a typical phase diagram, like that depicted in Fig.l, is determined ex- 
perimentally by preparing an array of samples of varying composition at room 
temperature and then heating each to temperatures above that required for com- 
plete homogeneous melting. Upon cooling, detection of a first order transition 
is possible by simple calorimetric methods. Much more sophisticated methods 
(including calorimetric, X-ray diffraction, etc.) must be used for detection of 
the second order transitions. Once the phase diagram is complete, the fraction 
of each component in a given phase can be extracted directly from the diagram 
at a specified temperature T and composition Co- 

For the sake of illustration, we begin with the simplest typical case of a 
two component (A and B) system existing in two phases, say a and (3, as 
depicted in Fig. 4. If such a system is kept under constant pressure, depending 
on temperature, we expect it to be either in the a or the (3 phases or in both. 
It is of interest to know the fraction of a given phase in the overall system, but 
it is not always possible to measure this directly. Normally, only the overall 
composition Cq is known. Although for a single-phase region on the phase 
diagram the composition is Co, for a two-phase region, the composition of each 
phase is obtained by first placing a horizontal line I through the point (Co,T) on 
the diagram. At each end, this horizontal line intersects the boundary separating 
the two phase region from a single phase region. The two-phase region is always 
bounded by two single-phase regions. Any larger diagram, like that depicted in 
Fig.l, is assembled of alternating fragments of this kind. The composition of a 
given phase is then obtained as an abscissa of the above mentioned intersection 
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A & Co % B 

Companion tyrfru B) 



Figure 4: The lever rule construction 

point. Thus, the compositions of the a and (3 phases are C a and Cp, respectively 
as depicted in Fig. 4. 

We would like to connect this information with the Landau order parameter 
ip. To this purpose, we introduce the following notation. Let Ip = \Cp — Co| 
and l a = \Cq — C a \ so that Ip + l a = \Cp — C a \ = I. With this notation, we 
let the total fraction w a of the a phase be w a — -j and accordingly, wp — l -f. 

By construction, w a + wp = 1 so that the Landau order parameter ip can now 
be defined as ip — \w a — wp\. To check that such a definition makes sense, 
it is sufficient to notice that since at criticality Cp — C a — Co, we obtain 
w a ~ wp = 1/2 in accord with Section 1. Such a picture assumes existence of a 
well defined single critical point as discussed in Rcf. [6], page 257. On the next 
page of this reference one finds the following statement: "Strictly speaking, 
there can be said to be two phases only when they exist simultaneously and 
in contact, that is at points lying on the equilibrium curve. It is clear that 
the critical point can exist only for phases such that the difference between 
them is purely quantitative, for example a liquid and a gas ..." Applying this 
statement to Fig. 4, one can say that two phases (co)exist as long as I ^= 
and the critical point corresponds exactly to the case when I — 0. Next, on 
the same page we read "...such phases as liquid and solid (crystal), or different 
crystal modifications of substance, are qualitatively different, since they have 
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different internal symmetry. It is clear that we can say only that a particular 
symmetry... element exists or does not exist; it can appear or disappear only as 
a whole, not gradually. In each state the body will have one symmetry or the 
other, and so we can always say to which of the two phases it belongs. Therefore 
the critical point cannot exist for such phases, and the equilibrium curve must 
either go to infinity or terminate by intersecting the equilibrium curves of other 
phases. " 

The background information we have just supplied is sufficient for under- 
standing of the CuZn phase diagram depicted in Fig.l. We begin with the re- 
gion where Zn concentrations are less than 31.9% and temperatures T > 900°C 
. In this domain the horizontal (temperature) line intersects two coexistence 
curves so that locally the picture looks exactly the same as in Fig. 4. Under 
such conditions we have a phase coexistence between the CuZn liquid and the 
Cu Zn solid alloy, whose crystal structure a is that for Cu, i.e. fee, Fig. 2a). 
By lowering the temperature below 900°C and increasing the concentration of 
Zn we observe the phase coexistence between two solid phases a and (3 where, 
according to the Table 1, the (3 phase is bec, Fig. 3a). This fact is in accord 
with results of Table 3 indicating that the free energy of the bec lattice is lower 
than that for fee so that such a structure can exist only at lower temperatures 
in accord with empirical observations discussed in the Introduction. Next, by 
going to still lower temperatures we obtain the (3 phase in coexistence with the 
a phase. The (3 phase is made of interpenetrating cubic lattices as depicted in 
Fig. 3. Such observation, again, is in accord with Table 3 (which tells us that 
the cubic lattice has still lower free energy than the bec lattice). The dashed 
line in this range of concentrations and temperatures represents the line of the 
second order phase transitions in accord with quotations from Landau and Lif- 
shitz book, Rcf.[6], stated above 35 . Exactly the same arguments are applicable 
to the next, (3 + L, (3 — and {(3 + j) — (0 + 7') portions of the phase diagram so 
that the dashed line (extending to concentrations just below 60 % according to 
Ref. [6] ) is still a line of second order phase transitions where the Lifshitz theory 
discussed in the Introduction is applicable. Notice that if for such concentra- 
tions of Zn we raise the temperature, then we would reach the curved (3 triangle 
which looks exactly like the curved a triangle but is located strictly below in 
the temperature range. This is again in complete accord with results of Table 
3. Next, notice that the 7 phase is made essentially of the body-centered cubic- 
type lattice and in this range of concentrations and temperatures above 700°C 
the situation is analogous to that for the earlier discussed a phase. Below 700°C 
and for concentrations above 68% the diagram exhibits apparent complications 
as depicted in Fig. 5. 

35 The lever rule depicted in Fig. 4 can be extended to cover this case by imagining the 
dashed line in Fig.l opening up a little bit initially, thus forming something like an eye, and 
then, finally closing up. In this case we always would get Cp = C a = Co and w a = wp = 1/2 
in accord with Landau and Lifshitz [6],Chr.l4. Under such conditions the order parameter ip 
is either or 1 in accord with above cited quotations from Ref. [6]. This causes no problems, 
however, in view of Eq.(1.3). 
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Figure 5: Fragment of Fig.l exibiting the most complicated phase behavior 



Although description of these complications lie beyond the scope of Lifshitz 
theory, we would like to discuss them briefly for the sake of completeness and 
in view of likely developments in the future. 

Since the rest of the phase diagram, Fig.l (in the range of concentrations 
above 80% and temperatures below 600°C), is simpler, we would like to discuss 
it first. Provided that we make formal substitutions of the type: a ^± e and 
/3 ^ r], this portion of the diagram looks exactly the same as the previously 
discussed portion including a + L, a, a + (3, (3 + L and (3 phases. This fact 
should not come as a total surprise in view of the fact that Cu and Zn have 
different types of lattices. According to the Table 1, they are of fee and hep type 
respectively. Nevertheless, the phase diagram does contain a surprise. It comes 
from the apparent contradiction between the results of our Table 3 requiring the 
free energies of fee and hep type crystals to be the same and obvious asymmetry 
of the CuZn phase diagram with respect to temperature axis. Earlier, we noticed 
that, indeed, the results of computer simulations support the analytical results 
of our Table 3. To resolve the apparent paradox we have to take into account 
several additional facts. First, we have to take into account that computer 
simulations as well as our calculations summarized in Table 3 have been made 
with account of the fact that the packing fractions (Appendix C) of both hep 
and fee lattices are the same. This is obviously correct for lattices whose vertices 
contain identical atoms modelled by the hard spheres and arranged in such a 
way that these spheres touch each other. In reality, however, metals are made 
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of atoms whose packing behavior may or may not be well modelled by the 
hard spheres. In particular, the literature data on pure metals, Ref. [88], page 
260, indicate that only Li, Na and Sr may form standard (Fig. 2, c/a ~ 1.633) 
hep lattices. Even for these metals such lattices are not the most common: 
for Li and Na the most common form is bec while for Sr is fee. In our case 
Cu has a stable fee lattice while Zn has the elongated (Fig. 2, c/a ~ 1.86) fee 
lattice. The fact that Li and Na exist in two modifications: bec and hep and 
that, experimentally, the bec lattice is more stable than hep, can be used for 
explanation of the portion of the CuZn phase diagram where the transition 
from 7 + L to 7 + 5 phase is depicted since 7 is of bec-type lattice while S is of 
hep-type. This fact is in apparent contradiction with the order of appearance 
of phases of lower symmetry with decreasing temperatures as discussed in the 
Introduction. Such an order is characteristic only for the phase transitions of the 
second order however (in accord with Ref. [6]) while a horizontal line between 
the 7 + L and 7 + S phases is characteristic for transitions of the first order 
as demonstrated experimentally in [7]. Moreover, in comparing our results 
against experimental data one should keep in mind that the results of Table 
3 are for the combination [3T while thermodynamically we have to look at T 
at different temperatures. When it comes to hard spheres, it is clear that the 
equality P x T fcc = P 2 T hcp holds if (3 1 = j3 2 implying T fcc = T hcp in accord 
with computer simulations [86]. However, for metals such as Cu and Zn, in 
view of what was said about their crystal structures, it is sufficient to require 
1 ^ cc (Cu) = j3 2 !F hcp {Zn) = N with Af being some constant. This leads us to 
the result: 

T fcc {Cu) = ^-^ hc P(Zn). 

Since the data presented in Fig.l imply that T\/T 2 > 1 this implies that 
T fcc {Cu) > T hcp {Zn) . This fact formally explains the source of the ap- 
parent temperature asymmetry of the phase diagram, Fig.l. Clearly, such an 
explanation is purely formal and should be replaced by quantum mechanical 
calculations. This task, however, is beyond the scope of this work. 

In looking at the phase diagram, Fig.l one should keep in mind how such 
diagram was obtained experimentally (we mentioned this already). That is, 
experimentally one cannot move horizontally across the diagram. This means 
that the respective free energies also should be compared with each other only 
vertically, i.e. temperature changes for the fixed composition. 

Now we are ready to provide our final comments regarding Fig. 5, i.e. we 
would like to discuss qualitatively the region of concentrations above 60% but 
below 80% and temperatures below 700°C. This region contains 3 triple points: 
two of them, P, at temperature slightly below 700°C, and the third at temper- 
ature slightly below 600°C. These tripple points are called peritectics. They 
are interesting because in both cases we have situation when either two phases 
are transformed into one (7 + L ^ S and S + L ^ e) or, as in the case of point 
E (called the eutectic), one phase 8 is transformed into two phases 7 and s. 
Although the Landau theory can be extended, in principle, to describe (locally) 
situations depicted in Fig. 5, attempts to describe the entire phase diagram, 
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Fig.l, quantitatively using the same model [12] thus far have not produced the 
desired results to our knowledge. 



Appendix A. Some facts about lattices, especially 3 dimensional. 

A lattice A is a free abelian group given by 

A = Zvi + ... + Zv n = {aivi + ... + a n v n : Oj e Z, i = 1, n}. (A.l) 

where Vi,...,v n is a basis in Euclidean space E™. Hence each lattice is the 
universal covering space for some n— dimensional torus T n . The fundamental 
parallelotope P(A) of the lattice A is given by 

P(A) = {Ai«i + ... + X n v n : < Xi < l,i = 1, ...,n}. (A.2) 

In solid state physics the fundamental parallelotope is known as primitive cell 
[13]. If the crystal lattice is composed of atoms of several different types, it 
is required that the atoms of a given type under basic translations are sent to 
atoms of the same type. Hence a lattice as a whole is a union of primitive lattices 
(e.g. see Fig. 2) so that the symmetry of the crystal lattice as a whole does not in 
general coincide with that for a specific primitive lattice. To distinguish between 
primitive lattices the concept of a volume for P(A) is useful. It is defined by 

volP(A) = \det(v!,...,v n )\ . (A.3) 

If the vector basis is changed by the unimodular transformation, volP(A) re- 
mains unchanged. In solid state physics such equivalent lattices (primitive cells) 
are called Bravis lattices. Since not all Bravis lattices are connected by the uni- 
modular transformation, in 3 dimensions there are 7 different crystal systems 
and 14 Bravis lattices associated with them. These lattices differ from each 
other by the length ratios between their edges and by the angles between these 
edges. The Bravis lattice of any kind has symbol P (for primitive) in solid state 
physics. The body-centered and the face-centered lattices are not primitive and 
are denoted by I and F respectively. Let G n (A) denote a (crystallographic) 
group of isometries of E" which maps A to itself. Fedorov, Schoenflies and Bar- 
low [89] have independently established that for n = 3 there are 230 such groups 
of isometries. Bieberbach has looked at the same problem in n-dimensions [90]. 
The 3rd column of Table 1 provides information about the particular catalog 
number (#) for the space group G${A). Every group G${K) is a composition 
of some orthogonal transformation a and translation T, i.e. Vg € C?3(A). We 
should have g = T o a. The orthogonal transformation is called rotation if 
deta = 1 and rotatory reflection if det a = —1. In general T o a ^ a o T. 
Clearly, a belongs to a point group H^(A) in the sense that for every Bravis 
lattice Ab it is a stabilizer, i.e. a As = Ab- As a group, H^(A) is finite and 
each of its finitely many elements is a rotation about some axis going trough 
the (arbitrarily chosen, e.g. coinciding with one of the lattice sites) origin O 



59 



by the angle which is an integer multiple of either tt/3 or 7r/2 36 . Let C n be a 
rotation which corresponds to the angle 2ir/n, then a rotatory reflection is given 
by S n = a o C n , i.e. it is a superposition of a rotation followed by reflection 
a (rotation by the tt angle) in the plane perpendicular to the axis of rotation. 
In particular, the inversion i — a o C2, while the reflection itself is given by 
aoE — E 7 where E is the unit element of the cyclic group C. The simplest group 
after the cyclic group C n is the group D n consisting of all rotations transforming 
regular n— sided prism into itself. This group has one n— th order C n axis and 
n second order axes perpendicular to it thus containing altogether 2n elements 
[9]. Next, we need the group T which consists of all rotations transforming a 
tetrahedron to itself. Analogously, the group O (the octahedral group) contains 
all rotations leaving the cube invariant. To bring all these facts to the level at 
which the symbolics of Table 1 can be explained, we need yet a couple more 
complications. First, we introduce a group is a symmetry group of a 

regular n— gonal pyramid containing one n— th order axis C n coinciding with 
the height of the pyramid and n vertical reflection planes that pass through C n 
axis. Reflections in these planes are elements of C nv . It can be shown [9] that 
this group is isomorphic to D n . Second, the group D n h is the symmetry group 
of a regular n— gonal prism differs from earlier introduced D n by the presence of 
additional n vertical reflection planes (just like C nv ) causing extra reflections in 
these planes. Third, the group T h is related to T cls Sji is related to C„. Lastly, 
the group O n is related to O in the same way as Th to T . 

Since each of 230 lattice space groups is made the same way as some compo- 
sition of the point group symmetry followed by translational symmetry, one can 
arrange these groups either in ascending or descending level of symmetry. In 
particular, to explain the data in Table 1 we need only lattices with point groups 
associated with cubic and hexagonal symmetries. The needed information can 
be arranged for instance in the following way 



The Bravis lattice 


The point symmetry 


The crystallographic symbol 


Cubic ( P, I ) 


O h ,T d ,0,T h ,T 


m3m, 43m, 432, m3, 23 


Hexagonal (P ) 




6/mmm, 6m2, 6mm, 622, 6/m, 6, 6 



Here the elements of point group symmetry are arranged from left to right 
in descending order (with the leftmost being the most symmetric). The crystal- 
lographic nomenclature symbols follow the ordering of point symmetry groups. 

All that was said about the direct lattice can be said also about its dual. To 
this purpose let x = {x\, ■ ■■,x n ) and y ~ (y\, ...,?/„) be two vectors of E"such 
that their scalar product x ■ y = x Y y\ + ... + x n y n .Then, the dual (or reciprocal) 
lattice A* in E™ is defined by 

A* = {x e E n : x ■ y e Z Vy e A}. (A.4) 

36 This can be easily understood using two dimensional plane R, 2 as an example. In Section 
5 we have mentioned that in case of Euclidean geometry such plane can be tesselated either 
by squares or by triangles. The three dimensional space can then be viewed as R 2 X R. 
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disc(A) = (2n) n — . From the definition of the matrix Aij it follows that 



The dual basis v\,...,v^ is defined in such way that v* ■ Vj — % so that the dual 
lattice can also be presented as 

A* = Zv\ + ... + Zv* n (A.5) 

In solid state physics [13] instead of the requirement x ■ y € Z in Eq.(A4) the 
alternative requirement x ■ y = 2irZ is used. It is motivated by the Fourier 
analysis and the Bloch theorem [11]. It can be easily demonstrated that vol (A) • 
vol(A*) = 1 while in the solid state physics one gets accordingly: vol (A) ■ 
vol(A*) = (2n) n . A lattice A is called integral if A C A* . This means that 
Va:, y € A one has x ■ y € Z . If v\, v n is the basis of such integral lattice, then 
the matrix = Vi ■ Vj is an integral matrix. In this case the discriminant of 
lattice A is defined by 

disc(A) — dct(A). (A.6) 

But [vol(A)] 2 — dot A and, therefore, [vol(A)] 2 — disc(A) = j~ \ . Again this 

' ' vol(A*) 

result should be amended in the case of solid state physics where we have instead 

vol(A) 

vol{A*) ' 

each vector Vi has its own decomposition : Vi = (vn, Vi„). Accordingly, for 
the integral lattice any vector x of such a lattice is made of matrix products of 
the type: x = £M with £ being a vector with integer coefficients: £ = ...,£„) 
and M being a matrix: M — — {vij}. Let M T be a transposed matrix. 
Then, the earlier defined matrix A (called the Gram matrix [52]) is made of M 
and M T according to the rule: 

A = MM T (A.7) 
Based on this definition, the norm N(x) of the vector x can be defined as: 

N A (x) = x ■ x = £A£. (A.8) 



Appendix B. Number-theoretic aspects of lattices. 

Based on the information presented in the Appendix A, we consider gener- 
alization of CM, discussed in Section 5, to multidimensional lattices. To this 
purpose consider the following eigenvalue equation 

xvi = auV! + ... + a ln v n (B.l) 
xv 2 = a 2 iv x + ... + a 2n v n 



xv n = a n ivi + ... + a nn v n 
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where the matrix elements G Z. This equation can be looked upon from 
different directions. For instance, if we consider it as an eigenvalue equation, 
i.e. 

det(a:I - A) = (B.2) 
then, written explicitly, we obtain an n-th order algebraic equation 

x n + a n _ 1 x n - 1 + ... + oo = (B.3) 

with integer coefficients. By definition, n solutions {xi} of such equation arc al- 
gebraic integers [91,92]. For instance, ±1, ±i ( the Gaussian integers discussed 
in Section 5) are solutions of the equation x 4 = 1 while j = |(— 1 + \f~~ 3) (also 
discussed in Section 5) is one of the solutions of the equation a; 3 = 1. In both 
cases solutions are roots of unity C{ n )k = cx p(±ik^-), k = 0, ...,n — l(for n = 4 
and 3 respectively). Solutions C( n )k by construction are independent of each 
other and can be considered as basis vectors (instead of Vi, v n ) in some vector 
space. This analogy can be made more precise. Let Xi be one of the solutions 
of Eq.(B.3). Substitute it into Eq.(B.3) and rewrite the result as 

-ao = Xidi + ... + x"a„ (B.4) 

(perhaps, with a n = 1). Then, the analogy with Eq.(A.l) becomes complete 
if we identify (x^, xf) with (vn, ...,v in ). From this analogy, the norm, the 
trace and the discriminant can be defined as well. There are some differences 
though connected with the Galois group symmetry. Since these topics are 
not immediately connected with results of the main text, the interested reader 
can look up these things in number theory literature. What is important for 
us, however, is the fact that only for quadratic number fields does the norm, 
Eq.(A.8), coincide with that defined in number theory. For higher order fields 
this is no longer true [92]. Therefore, two dimensional results presented in 
Section 5 are not immediately generalizable to higher dimensions. 

Fortunately, there is an alternative approach bypassing the difficulty just 
described, which we would like to describe now. We begin with defining the 
Epstein zeta function Z(A, s) is defined as 

Z{A,s) = \ Y, [Na(x)}-\ (B.5) 
?eZ"\o 

where Na(x) is defined by Eq.(A.8). 
Using the known identity 

oo 

x- s T(s, x) = j dyy 3 - 1 exp(-xy) (B.6) 
o 

we can rewrite Eq.(B.5) as 

oo 

A(A,s) = \J dyy'-\(B{A;y) - 1) (B.7) 
o 
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where A(A, s) = ir s T(s)Z(A, s). The theta function 

®(A;y) = exp(-nN A (x)y) (B.8) 

can be represented alternatively as follows 

oo 

&(A;y) = J2<n)q n , (B.9) 

n=0 

where q — exp(— irny) and a(n) is the number of integer solutions of the equation 
Na(x) = n. These results are used in Section 6. 



Appendix C. Exterema of the positive definite quadratic forms and 
the problem of closest packing of spheres. 

Although methods discussed in Section 5 are capable of providing (indirectly) 
the information about the best (closest) possible packing of hard discs, they are 
not immediately generalizable to higher dimensions. Because of this, we sketch 
here an alternative approach. To this purpose, let us consider a positive definite 
quadratic form 

f(Sl, &) = <*ll£? + 2012^2 + 022^ = HA£ (C.l) 

whose determinant D = a\\aii — a\ 2 ■ Consider now a subset 1Z of £— plane 
such that £2) < K with K being some positive number. Following Cassels, 
Ref.[93], we would like to demonstrate that, provided that, K > (4£)/3)5, the 
subset 1Z contains some point (£1,^2)) other than the origin (0,0), with integer 
coordinates 37 . This result can be understood based on the following chain of ar- 
guments. First, let us assume that inf /(mi,u 2 ) = an with U\,U2 being integers 
not both equal to 0. Clearly, every quadratic form can be rewritten as 

/(£i,£ 2 ) = M£i + — £ 2 ) 2 + — & (C2) 

an an 

so that by construction, f(ui,U2) > an. Second, without loss of generality, we 
can demand that f{u\, 1) > an. If, in addition, we assume that 
then, we obtain the following chain of inequalities 



1 ail 



< 



an</(wi,l)<^ i + — (C.3) 
4 an 

leading to 

a 2 u < ^ (C4) 



37 In d-dimensions the proof of this fact is known as the Minkowski theorem. It plays a 
major role in geometric number theory and the theory of polytopes and oriented matroids. 
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and, finally, to 

inf/(^^ 2 )<A/^. (C.5) 




Consider a special case of the quadratic form £ 2 ) — m (£1 + £1^2 + £2) • ^ 
has determinant D = ^j- ■ So that inequality (C.5) is reduced to the identity 
: m — m. For any other quadratic form the sign < in (C.5) should be replaced 
by < [94]. The quadratic form ^ + + £2 corresponds to the hexagonal 
lattice [52,94] whose matrix M is given by 



This matrix gives rise to the matrix A = MM T . Two of such matrices are 
equivalent if they are connected by the unimodular transformation made of 
matrix U with integer entries whose determinant is one, i.e. A' = UAU T . This 
observation permits us to use as well the matrix [52], page 43, 

M =(4^) <C7) 

made of two basis vectors (ideals) as discussed in Section 5. 

Consider the obtained results from a slightly different angle. Let (in any 
dimension) 

a = inf£A£ (C.8) 

for £ e Z"\0 . Let each lattice point be filled with a sphere of radius r and 
let the radius r be such that spheres at nearby points touch each other. Then, 
2r < ^fa. In view of Eq.(A.6), the volume of the fundamental parallelotope, 
vol(A) = V dct A. If the volume V n of the n dimensional sphere is given by 

- r(i+f ) r ™ = a nr n , then the most optimal packing fraction q n can be 
defined as 

where the constant 7,? is defined by the following inequality 

a< 7n (detA)". (CIO) 

According to Eq.(C.5) we obtained in two dimensions : 7 2 = -j=- This produces 
at once q 2 = ^ Jj- = ^= = -J= = 0.9069 in accord with Ref.[52], page 110. In 
the case of 3 dimensions, the analog of the hexagonal two dimensional lattice is 
the face centered cubic lattice [52] (fee) whose quadratic form is known to be 

/ (£1, £ 2 , £3) - i\ + il + il + ?i? 2 + £ 2 ? 3 + £ 3 £i- (en) 

Calculations similar to those leading to inequality (C.5) were made originally 
by Gauss and can be found in Ref. [61,93,94]. They now lead to the following 
inequality 

a<(2dctA)* (C.12) 
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thus yielding 7 3 = 2i and q 3 = ^-^\f2 = ^= = 0.7405. This value is not 
unique, however, since the hexagonal closed packed (hep) lattice also possesses 
the same value for q 3 . 



Appendix D. Representative calculation of determinant without zeta 
function regularization 

The quantum harmonic oscillator is a benchmark example of quantum me- 
chanical calculations in any course of quantum mechanics. Its eigenvalue spec- 
trum is known to be: E n = Huj(n+ i), where u> is the classical frequency of the 
oscillator and h is Planck's constant. Given this result, the partition function 
Z(f3) at temperature /3 _1 is obtained in a standard way as 

Using this result the free energy T is obtained as 

0F=- lnZ(/3) = (3— + ln(l - exp(-/?fiw))- (D.2) 

We would like to reproduce this simple result now using Feynman's path integral 
approach to quantum and statistical mechanics. To this purpose, following 
Feynman [95], we write for the free energy the following path integral (H = 1) 



cxp(-/3^)= J D[q{r)] C M-\J dT[q 2 +0J 2 q 2 ]) (D.3) 
<z(0)=g(/3) 

where q = J^(t). Since for given boundary conditions for g(r) we can write 

2 

drq 2 = - j drq(^)q (D.4) 



o 



this fact allows us to rewrite the path integral, Eq.(D.3), in the following equiv- 
alent form 

l l 

exp(-#F)= J D[q(T)]exp(~ J dr J dr' qA(T,r')q) (D.5) 

q(Q)=q(l) 

where the operator A(t,t') = Jpr + ll> 2 /3 2 ^J S(t — r'). By expanding q(r) in 
Fourier series with basis made of eigenfunctions of such operator the calculation 
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of the path integral is reduced to the calculation of the Gaussian path integrals 
of the standard type 

oo 

f dxcxp(-^x 2 ) = (D.6) 

J ^ V a n 

— oo 

where a n = 47r 2 n 2 + tj 2 /3 2 ,n = 0, ±1, ±2, ... Thus, we obtain 

oo 

detA=J|a„ (D.7) 

— oo 

and, therefore, the free energy is formally given by 

f3T =ilndct A + const (D.8) 

where the actual value of const is unimportant to us since it is temperature- 
independent and can be dropped since the free energy T is always defined with 
respect to some reference. To calculate the free energy we follow the method 
described in our earlier work [82]. To this purpose, we let lo 2 (3 2 — x 2 and 
consider 

d 1 °° 1 

^ lndet ^ + 2 £«T^ (D ' 9) 

n=l 

Then, formally we obtain 



In det A = I dx 



oj 2 /3 

— 2V ' 



n=l 



(D.10) 



where the lower limit of integration will be carefully chosen. By noticing that 

1 2x °° 1 
COth7TX= — + — (D.ll) 

TTX TT ' n Z + X* 
n=l 

we obtain 



d 1 x 

— In det ,4 = — coth- (D-12) 

dx 2 2x 2 v ' 



and, accordingly, 

W3 



, det A(u)B) 

In ., \ 

det A(0) 



= y*dxcoth| =w/3 + 21n(l-exp(-/3w)) (D.13) 



Finally, taking into account Eq.(D.8) we obtain (H = 1) 



$T + ln(l - exp(-/3w)) (D.14) 



in accord with Eq.(D.2). 
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